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ABSTRACT 

The international scientific community is actively 
engaged in assessing ice sheet and alpine glacier 
fluctuations at a variety of scales. The availability 
of stereoscopic, multitemporal, and multispectral 
satellite imagery from the optical wavelength 
regions of the electromagnetic spectrum has greatly 
increased our ability to assess glaciological con- 
ditions and map the cryosphere. There are, how- 
ever, important issues and limitations associated 
with accurate satellite information extraction and 
mapping, as well as new opportunities for assess- 
ment and mapping that are all rooted in under- 
standing the fundamentals of the radiation 
transfer cascade. We address the primary radiation 
transfer components, relate them to glacier 
dynamics and mapping, and summarize the analy- 
tical approaches that permit transformation of 
spectral variation into thematic and quantitative 
parameters. We also discuss the integration of 
satellite-derived information into numerical model- 
ing approaches to facilitate understandings of 
glacier dynamics and causal mechanisms. 


2.1 INTRODUCTION 

Concerns over greenhouse gas forcing and globally 
increasing temperatures have spurred research into 


climate forcing and the responses of various Earth 
systems. Complex geodynamics regulate feedback 
mechanisms that couple atmospheric, geologic, 
and hydrologic processes (Molnar and England 
1990, Ruddiman 1997, Bush 2001, Zeitler et al. 
2001, Bishop et al. 2002). A significant component 
in the coupling of Earth’s systems involves the 
cryosphere, as glaciological processes affect atmo- 
spheric, hydrospheric, and lithospheric responses 
(Bush 2000, Shroder and Bishop 2000, Meier and 
Wahr 2002). Scientists have thus recognized the 
significance of understanding glacier fluctuations 
and their use as direct and indirect indicators of 
climate change (Kotlyakov et al. 1991, Seltzer 
1993, Haeberli and Beniston 1998, Maisch 2000). 

The international scientific community is actively 
engaged in assessing ice sheet and glacier fluctua- 
tions at the global scale (Zwally et al. 2002, Bishop 
et al. 2004a, Dyurgerov and Meier 2004, Kargel et 
al. 2005). Glaciological studies using remote sensing 
and geographic information system (G1S) studies 
indicate that recent glacier retreat and wastage 
appear to be global in nature (Dyurgerov and Meier 
2004, Paul et al. 2004, Barry 2006, Zemp et al. 
2006). It is essential that we identify and character- 
ize those regions that are changing most rapidly 
and having the most significant impact on sea 
level, water resources, economics, and geopolitics 
(Haeberli 1998). Many regions have not been ade- 
quately studied, however, and progress in scientific 
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understanding requires detailed analysis of param- 
eters such as the number of ice masses, land ice 
distribution, changes in glacier length and volume, 
mass balance, dynamic response time, and ice flow 
dynamics. 

Haeberli (2004) and Zemp et al. (see Chapter 1 of 
this book) indicated that a multitiered approach to 
this problem is required, involving measurements 
across environmental gradients, mass balance stu- 
dies, regional assessment within different climatic 
zones, long-term observations of glacier length 
changes, and systematic multitemporal land ice 
inventories. This approach integrates field-based 
glaciological research, numerical modeling, and 
remote sensing and GIS-based investigations. 

Careful in situ mass balance and other ground 
reference studies are an essential component of 
scientific investigations of glacier responses to 
climate change. However, McClung and Armstrong 
(1993) suggested that even detailed studies of well- 
monitored glaciers do not permit characterization 
of regional or global mass balance trends. New 
developments in remote sensing are required to 
bridge observations at a variety of scales, from 
the individual glacier to the global picture (Haeberli 
1998, Bishop et al. 2004a, Dyurgerov and Meier 
2004, Raup et al. 2007). 

Glaciological research has increasingly used data 
from microwave, optical, and light detection and 
ranging (LiDAR) sensors aboard various space- 
based and airborne platforms. Collectively these 
data and new methodologies have been used to 
map glacier boundaries, ice flow velocity fields, 
and estimate other glaciological parameters such 
as the equilibrium line altitude (ELA) and mass 
balance. The availability of stereoscopic, multi- 
temporal, and multispectral optical satellite 
imagery has greatly increased our ability to assess 
glaciers and map their distribution (Bishop et al. 
2004a). There are important opportunities and 
limitations associated with accurate satellite 
information extraction and mapping that are 
rooted in the fundamentals of the radiation transfer 
cascade. 

The purpose of this chapter is to examine the 
fundamentals of remote-sensing science that serve 
as a basis for information extraction and glacier 
assessment. We address the primary components 
of radiation transfer theory and relate these to 
glacier dynamics and mapping. Finally, the use of 
satellite-derived information in numerical modeling 
is presented for understanding causal mechanisms 
responsible for ice mass fluctuations. A list of 


parameters and their definitions can be found 
toward the end of this chapter (Section 2.9, p. 47). 

2.2 RADIATION TRANSFER CASCADE 

Characterization and mapping of glaciers using 
multispectral and multitemporal satellite imagery 
pose a variety of scientific and technical challenges 
(Bamber and Kwok 2004, Bishop et al. 2004a). 
Numerous physical parameters of the atmosphere 
and the landscape directly and indirectly govern 
radiation transfer (RT) processes. Radiance meas- 
ured at the sensor, therefore, does not accurately 
represent the reflectance or emission of radiance 
from the landscape surface. Image spectral varia- 
tion may be due to variations in atmospheric prop- 
erties or topography, rather than surface material 
properties. The chain of RT processes, which 
integrates the influence of multiscale atmospheric, 
topographic, and surface matter-energy interac- 
tions, is complex. 

For glaciologists interested in using satellite 
imagery, knowledge of the RT components is neces- 
sary because their influence governs spectral 
variability in imagery. Spectral variability due to 
coupled atmosphere-topography-surface interac- 
tions is problematic for applications that rely on 
calibrated surface radiance values. Sensor system 
characteristics also affect spectral variation in 
imagery (e.g., spatial, spectral, and radiometric 
resolutions). Therefore, Earth scientists must make 
decisions regarding radiometric calibration, image 
enhancement, and choice of pattern recognition 
algorithms. Depending upon environmental con- 
ditions. RT components may or may not signifi- 
cantly influence image spectral variability. 
Unfortunately, the influence of individual RT com- 
ponents (and the entire cascade) on glacier assess- 
ment and mapping has not been adequately 
researched. 

2.2.1 Solar irradiance 

The RT cascade begins with the energy output of 
the Sun and ultimately traces through to a sensor 
aboard a mobile platform. Numerous processes in- 
fluence the magnitudes of emission, absorption, 
transmittance, and reflection of light coming from 
the Sun. Besides the obvious and important control 
of the magnitude of solar irradiance by latitude, 
time of day (solar elevation), and season (annual 
orbit and current obliquity), solar irradiance also 
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varies as functions of: (1) solar dynamics (magnetic, 
sunspot, and solar flare cycles related to the solar 
dynamo and convection); (2) Earth’s orbital 
eccentricity and rotational dynamics (including 
the changing orbital ellipticity and distance from 
the Sun to the Earth, rotational obliquity, and 
obliquity phase with respect to Earth’s orbital 
phase); and (3) changes in the radiation field on 
the Earth’s surface due to atmospheric scattering 
and absorption (including changes in infrared 
opacity due to greenhouse gas and volcanic emis- 
sions). Consequently, there are changes in irradi- 
ance that are wavelength dependent ((1) and (3) 
above), and others that are wavelength independent 
((2) above); and changes that are strictly cyclic and 
modulated at many frequencies ((1) and (2) above), 
and others that are not cyclic beyond the annual 
primary cycle but have either secular trends 
(anthropogenic greenhouse gases) or episodicity 
(natural greenhouse gases and anti-greenhouse vol- 
canic aerosols). Each of these is important from the 
perspective of long-term climate change and glacial 
conditions, but from the remote-sensing image 
analysis perspective these variations are considered 
relatively insignificant because variability is either 
of a long period (thus having small interannual 
variability) or has a small amplitude. 

The spectral exitance (M) of the Sun and Earth 
can be approximated using Planck’s blackbody 
equation: 

^® J A) = A ^e xpC2/(A 1 re) _ 10 ^ (2- 1 ) 

where T e is the temperature of the Sun, and C) and 
C 2 are the first and second radiation constants. The 
exoatmospheric spectral irradiance ( E ° ) at the top 
of the Earth's atmosphere can be approximated 
such that 



Figure 2.1. Comparison of exoatmospheric spectral 
irradiance from the composite SMARTS2 spectrum 
and a blackbody approximation. The SMARTS2 spec- 
trum is a composite spectrum based in part on the 
WRC85 spectrum and balloon and satellite observa- 
tions (after Gueymard 1995). 

2.2.2 Surface irradiance 

Radiation is transmitted through the atmosphere, 
and atmospheric constituents and topographic fac- 
tors dictate the magnitude of radiation that reaches 
the surface. Surface irradiance ( E ) is a composite of 
three downward irradiance components that are 
wavelength (A) dependent such that 

E(X) = E b (X)+E d (X) + E t (X), (2.3) 

where E h is direct beam irradiance, E d is diffuse 
skylight irradiance ( E d ), and E, represents adjacent 
terrain irradiance. Individually and collectively 
these irradiance components govern a variety of 
surface processes such as ablation, and also dictate 
the magnitude of surface reflectance and emission. 
Depending upon environmental conditions, these 
components have varying degrees of influence on 
image spectral variability. 


E°( A) 


M(T e ,X) A Q 

7 r 


( 2 . 2 ) 


where A 0 is the area of the solar disk and d es is the 
Earth-Sun distance which is a function of time. 
Variations in measured exoatmospheric spectral 
irradiance are caused by solar atmospheric interac- 
tions, and this can be compared with the blackbody 
approximation (Fig. 2.1). Standard exoatmospheric 
irradiance spectra can be used for the computation 
of spectral band irradiance values needed for the 
radiometric calibration of imagery. 


2.2.2. 1 Direct solar spectral irradiance 

Under cloudless skies the direct irradiance com- 
ponent constitutes the dominant term in equation 
(2.3). Much research has been conducted to model 
atmospheric transmittance functions accurately. 
The atmosphere attenuates direct irradiance pri- 
marily by gaseous absorption and molecular and 
aerosol scattering (Chavez 1996). These atmo- 
spheric processes are wavelength dependent and 
spatially and temporally controlled by changing 
atmospheric and landscape conditions. Total down- 
ward atmospheric transmission (T^) is a function of 
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Figure 2.2. Simulated atmospheric transmittance in the shortwave spectrum at the terminus of the Baltoro Glacier 
near K2 in Pakistan on August 14, 2004 at 10:00 am. Latitude: 35.692243° N; longitude: 76.165222° E; altitude: 
3,505 m; solar zenith angle: 36.19°; precipitable water: 1.0 cm; ozone: 0.35 cm. 


the total optical depth of the atmosphere which 
varies with solar zenith angle and altitude, and 
can be represented as 

T*(A) = T,(A)T fl (A)T 03 (A)V(A)T H2 o(A), (2.4) 

where T,. is Rayleigh transmittance, T„ is aerosol 
transmittance, Tq 3 is ozone transmittance, T gas is 
transmittance for miscellaneous well-mixed gases, 
and T HiG is water vapor transmittance. Atmo- 
spheric attenuation is highly variable with 
wavelength, with Rayleigh and aerosol scattering 
dominating at shorter wavelengths and water vapor 
absorption at longer wavelengths (Fig. 2.2). Direct 
irradiance into the surface is also governed by 
multiscale topographic effects. Local or microscale 
topographic variation is represented by the inci- 
dence angle of illumination (/) between the incom- 
ing beam and the vector normal to the ground, such 
that 

cos i = cos 8j cos 6, + sin 9 t sin 6, cos (</>, — </>,■), (2.5) 

where 9 t is the incident solar zenith angle, 4>j is the 
incident solar azimuth angle, 9 , is the slope angle of 
the terrain, and <j) t is the slope-aspect angle of the 
terrain. 

Estimation of i is possible with the use of a digital 
elevation model (DEM), and uncertainty in the 
estimate is related to the measurement scale, as 
subpixel-scale topographic variation is not 
accounted for. Values of cos i can be <0.0, indi- 
cating no direct irradiance due to the orientation of 
the topography. It is important to note that the 


incident solar geometry varies across the landscape, 
although this is usually assumed to be constant 
when working with individual image scenes (i.e., 
small-angle approximation). In addition, mesoscale 
topographic relief in the direction of </>,■ determines 
if a pixel is in shadow ( S ). This parameter value will 
be 0.0 or 1.0 depending upon the presence or 
absence of a cast shadow, respectively. Satellite 
imagery acquired in rugged terrain or with rela- 
tively large solar zenith angles will exhibit cast 
shadows. This can be accounted for by ray tracing, 
shadow detection, and shadow interpolation algo- 
rithms that alter cos i values appropriately (Dozier 
et al. 1981, Rossi et al. 1994, Giles 2001). Local 
topography and shadows cast increase the global 
spectral variance in satellite images, with a decrease 
in spectral variance within shadowed areas. Collec- 
tively, the aforementioned parameters define the 
direct irradiance component as 

E b ( A) = £'°(A)T l (A) cos iS. (2.6) 

In high-latitude and mountainous environments, 
this irradiance component can have a significant 
influence on glacier dynamics and image spectral 
variability. Local topographic variations in param- 
eters such as slope magnitude and azimuth lead to 
significant variation in E h due to ablation, supra- 
glacial fluvial action, and supraglacial debris trans- 
port processes. Furthermore, spatial variations in 
slope occur as a function of altitude (Fig. 2.3). 
Finally, the glacier altitude range determines the 
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Figure 2.3. Slope-altitude plots for four alpine glaciers in the Karakoram Himalaya. Geomorphometric analysis 
was conducted using a SRTM 90m digital elevation model. Average slope values were computed over a 50m 
altitude interval. Different magnitudes and altitude trends are depicted. Such variation increases when a smaller 
altitude interval is used. 


relative influence of atmospheric transmittance on 
this component. 

In general, this component is significant and 
highly variable across most mountain landscapes 
exhibiting moderate to extreme relief (Fig. 2.4). 
Such variability is present in satellite imagery and 
usually needs to be accounted for using topographic 
normalization techniques. Unfortunately, many 
methods of topographic normalization do not 
effectively remove this component from satellite 
imagery, as they do not inherently account for 
multiscale topographic effects. Automated algo- 
rithms for thematic mapping are sensitive to such 
spectral variations caused by E b and can classify 
steep slopes, cast shadows, moisture-laden debris, 
and surface water as the same feature. Conse- 
quently, classification accuracy can be highly vari- 
able depending upon the extent of topographic 
variability. Empirical spectral feature extraction 
(e.g., image ratioing and principal component 
analysis) techniques are frequently used to reduce 
these topographic effects, although they can still 
influence classification results depending upon what 
approach or algorithm is used. 

2.2.2.2 Diffuse skylight spectral irradiance 

Atmospheric scattering generates a hemispherical 
source of irradiance. This source can be simplis- 


tically represented as a composite including a 
Rayleigh-scattered component (E r ), an aerosol- 
scattered component ( E a ), and a ground- 
backscattered component (E g ) that represents inter- 
reflections between the landscape surface and the 
atmosphere, where 

E d (X) = E r (X)+E a (X)+E g (X). (2.7) 

Its accurate estimation is complicated by the fact 
that an anisotropic parameterization scheme is 
required. In general, diffuse skylight irradiance 
decreases with angular distance from the Sun. In 
addition, this irradiance component is also influ- 
enced by mesoscale hemispherical shielding of the 
topography. Consequently, only a solid angle of the 
sky will contribute to E d , and this angle will change 
as a function of pixel location and azimuth. In 
general, the solid angle will increase with altitude. 
It is frequently referred to as the sky view factor 
(Vf) in the remote-sensing and energy balance 
literature, and can be estimated using a DEM such 
that 

27T a / 

Vf = Y cos 2 6 max {<f>, d) — , (2.8) 

( j )—0 Z7r 

where 6 max is the maximum local horizon angle at 
a given azimuth, q !>, over a radial distance of d. In 
mountain environments exhibiting extreme relief 
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Figure 2.4. Simulations of local and mesoscale topographic influences on direct irradiance for the Mt. Everest 
region in Nepal. Grayscale values represent the collective magnitude of cos iS that is proportional to E b without the 
influence of the atmosphere. Simulation results were obtained using ASTER GDEM data and different solar 
geometry. A. Solar zenith and azimuth angles of 0 and 135°, respectively. B. Solar zenith and azimuth angles 
of 45 and 1 35°, respectively. C. Solar zenith and azimuth angles of 70 and 1 35°, respectively. D. Solar zenith and 
azimuth angles of 85 and 135°, respectively. 


(deep valleys), topographic shielding of skylight 
diffuse irradiance can be significant (Proy et al. 
1989). Consequently, the sky view factor is impor- 
tant in environments where climate variability and 
glacier erosion can generate moderate to extreme 
topographic relief (Fig. 2.5). It will be relatively 
insignificant in environments with minimal topo- 


graphic relief. Where significant, individual glaciers 
within a region receive variable amounts of energy 
from this component. Glacier fluctuations and mass 
movements alter the relief and ridgelines, such that 
glacier surfaces exhibit sky view variations with 
altitude (Fig. 2.6). Furthermore, difficulties inaccu- 
rately predicting the bidirectional reflectance distri- 
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Figure 2.5. Sky view factor over the Mt. Everest region in Nepal. Dark gray tones represent relatively high 
topographic shielding of diffuse skylight irradiance. Lighter tones signify more diffuse skylight irradiance reaching 
the landscape (isotropic assumption). These results were obtained using ASTER GDEM data. Examination of the 
spatial patterns clearly indicates that individual glaciers receive different amounts of energy and that altitude trends 
exist. 



Skyview Coefficient 

Figure 2.6. Sky view-altitude plots for four glaciers in the Karakoram Himalaya. Geomorphometric analysis was 
conducted using a SRTM 90 m digital elevation model. Average sky view coefficient values were computed over a 
50 m altitude interval. Different magnitudes and altitude trends are depicted. 
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bution function (BRDF) of land cover adjacent to a 
particular location generates uncertainty in the esti- 
mation of the ground-backscattered component. 

Such mesoscale land cover and topographic 
effects on this irradiance component are present 
in satellite imagery, although visual detection is 
difficult due to the dominance of direct irradiance 
and surface reflectance. In addition, lower spatial 
frequency components are difficult to detect and 
reduce, as topographic normalization approaches 
do not effectively address this irradiance com- 
ponent. 

More research is required to determine what the 
magnitude of this topographic effect is on image 
spectral variability. However, the magnitude can 
be deemed as potentially crucial in certain limited 
but important cases. The diffuse skylight irradiance 
component generally is of the order of 5-10% of 
direct irradiance, as any photographer knows, or as 
gauged from the ASTER VNIR-band DN values in 
typical mountain shadows compared with directly 
illuminated terrain. In heavily shadowed zones 
(such as in northern-aspect cirque glaciers in the 
Northern Hemisphere), the diffuse skylight com- 
ponent can approach the direct solar beam in its 
importance to the energy budget. The sky view 
factor in rough terrain, such as in Himalayan val- 
leys, can then account for blocking up to half, or 
even more, of the diffuse skylight. In general, rela- 
tive to direct beam irradiance, diffuse skylight and 
skyview shielding each become more important as 
relief increases (e.g., in narrow canyons) and as 
latitude increases (MacClune et al. 2003). 

2.2.2. 3 Adjacent terrain spectral irradiance 


where L u represents the radiance received at p\ 
from the luminance of p 2 ( L 2 ), 0\ and 0 2 are the 
angles between the normal to the terrain and the 
direct line of sight from/;! to p 2 , A p is the pixel area 
( p 2 ), and d represents the distance between p x and 
Pi- 

This equation can be used to estimate E, for any 
pixel by integrating over all of the pixels whose 
slopes are oriented towards a pixel of interest and 
where the line of sight is not blocked by the topog- 
raphy. High-altitude and extreme local relief areas 
can exhibit a strong adjacent terrain irradiance 
component due to highly reflective features such 
as snow and vegetation. 


2.2.3 Surface reflectance 


The magnitude of reflected and emitted radiance at 
the surface is determined by the conservation of 
energy, and is expressed by Kirchhoff ’s law as 

p(X) + a a (X) + T(A) = 1.0, (2.10) 


where, p, a„, and T represent reflectance, absorp- 
tion, and transmission, respectively. For opaque 
objects T = 0.0. Therefore, spectral reflectance 
can be approximated as 


PW = 


M(T S ,X ) 

E( A) 


( 2 . 11 ) 


where T s is surface temperature. Given a Lamber- 
tian surface, which reflects radiation equally in all 
directions, reflected surface spectral radiance (L) 
can be computed as follows: 

L(X) = p(X)(&j. (2.12) 


The irradiance components E h and E d interact with 
the terrain and land cover biophysical characteris- 
tics to generate an adjacent terrain irradiance com- 
ponent (E t ). 

This irradiance component is not generally con- 
sidered in remote-sensing and energy balance 
studies because it is assumed that its magnitude is 
relatively minor and it is a difficult parameter to 
estimate accurately. A first-order approximation 
was formulated by Proy et al. (1989) and assumes 
that reflected surface radiance is Lambertian. It is 
then possible to estimate the radiance received at 
any pixel by accounting for the geometry between 
two pixels ( p\ and p 2 ) such that 


L 12 = cos 9\ 


^ L 2 cos d 2 



(2.9) 


Satellite imagery can be used to estimate surface 
radiance and reflectance, although atmospheric 
correction is required. In general, this can be 
approximated as 

L( A) = (L°(X)-L p (X))/T\e r ,X), (2.13) 

where L° is at-satellite measured radiance, L p is 
additive path radiance from the atmosphere, and 
T 1 is upward total atmospheric transmittance that 
is a function of the reflected zenith angle (0,.) and 
wavelength. 

2.2. 3.1 Bidirectional reflectance distribution 
function and albedo 

Unfortunately, the Lambertian assumption is not 
valid on glaciers (see Chapter 3 of this book by 
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Furfaro et al.), and the anisotropic nature of surface 
reflectance must be characterized to permit accurate 
assessments of surface albedo and ablation. 
Numerous parameters influence reflectance aniso- 
tropy which is characterized by the bidirectional 
reflectance distribution function (f). It is defined 
as the ratio between reflected radiance and units of 
incoming surface irradiance. It is represented as 


0/ ■> 0 r , (p r , A) 


L(0i : (pi , 0 r ■ 0/-. A) 

E(O h cPn A) 


(2.14) 


where 0,. is the reflected azimuth angle. This func- 
tion characterizes bidirectional reflectance for all 
combinations of incident and viewing angles. 

In satellite remote sensing, a distinct, but closely 
related, parameter is the bidirectional reflectance 
factor ( R r ). This parameter is frequently utilized 
because it characterizes directional reflectance in a 
convenient unitless form. It is formally defined as 
the ratio between the reflected radiance in any par- 
ticular direction, and the radiance that would be 
reflected into the same direction by an ideal Lam- 
bertian surface, illuminated by the same source with 
the same incident geometry. It can be shown (Feng 
et al. 1993) that f. and R r are related by the simple 
relationship 

R,.(9 h 0,., <p„ A) = 7t f r (0 t , ^ 0 r , <p r , A). (2.15) 


This quantity is becoming a standard product for 
satellite-based radiometers capable of acquiring 
angular measurements of reflected radiation from 
the landscape over a scene (e.g., multiangle imaging 
spectroradiometer, MISR, and the polarization and 
directionality of Earth reflectance instrument, 
POLDER; see Diner et al. 1999). The values of 
R,. can be greater than one, for example, due to 
specular reflection. 

The aforementioned reflectance parameters are 
closely related to a parameter generally called 
surface albedo (a). It is a spectrally and hemispheri- 
cally integrated surface parameter that significantly 
governs surface energy balance. We distinguish 
between two types of albedo. Planar surface albedo 
( a psa ) is formally defined as follows: 

p2n /» 7 t /2 

^psa^ii / / VrfriQii fin &■> dQ r d(j) n 

Jo Jo 

(2.16) 

where = \n-/l r \, n is the normal vector, and 
= (6 r ,(p r ). It is the ratio between reflected and 
incoming radiation integrated over the hemisphere. 
The integration removes the dependence of albedo 


from the reflection geometry and therefore it is a 
global quantity. 

Spherical surface albedo (a ssa ) is obtained by 
averaging planar albedo over the irradiance hemi- 
sphere. It is defined as 

1 /■27T rn/ 2 

&ssa( A) —J I 0 J i®psa(O i> (pi- A) dO L d(pi, (2.17) 

where p,- = | n • O,- and O,- = (6,, 0 ( ). It is important 
to point out that for any surface in the cryosphere, 
Kirchhoff’s law of conservation during matter- 
energy exchange (assuming local thermodynamic 
equilibrium) dictates that emissivity and albedo 
are related by the simple relationship 
e s (A) = 1 — a ssa (X), respectively. 


2.2.4 Surface emission 


Upward longwave surface spectral radiance L 1 can 
be approximated as 


L t (A) 


£j(A) 


M(T S , A) 

7 r 


(2.18) 


Emissivity is a parameter that characterizes the 
efficiency with which the surface radiates energy 
compared with an ideal cavity or blackbody, for 
which e s = 1. 

If e s is known for at least one reference wave- 
length, as it is for snow and ice (Fig. 2.7), surface 
temperature can be computed as 

T s = AC 2 -' [ln( e ,(A)C 1 7r- 1 A- 5 Lf 1 (A) + l)]' 1 . (2.19) 


Since T s is invariant with A, it is then possible to 
calculate the emissivity spectrum if multispectral 
data are available. For most remote-sensing studies, 
however, other factors must be considered, compli- 
cating the inversion for T s (i.e., atmospheric 
effects). Consequently, at-satellite measured long- 
wave spectral radiance l/\ can be represented as 

L?(A) = TT(A)( e ,(A)M^) + (1- £ J (A))U l (A)7 r -'j 

+ Lj(A), (2.20) 

where T T is upward longwave atmospheric trans- 
missivity. Downwelling longwave spectral irradi- 
ance is L/ and l) p is upwelling longwave 
atmospheric path radiance. is reflected from 
the surface, which is usually taken to be a perfect 
diffuse (Lambertian) reflector, in which case reflec- 
tivity is simply given by (1 — £, 5 (A)) according to 
Kirchhoff’s law. 
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Figure 2.7. Thermal emissivity of water, snow, and 
ice. The uncertainties of emissivities estimated from 
ASTER are about ±1.5%. Spectral coverage does not 
extend to long enough wavelengths to exploit the 
spectral differences of snow and ice fully. The emissiv- 
ities calculated from reflectivity data measured by 
J.H. Salisbury at Johns Hopkins University (http:// 
speclib.jp/. nasa.gov/ scripts / lib /asp/ WaterResp.asp ) , 
using Kirchhoff's Law. The emissivity of natural rough 
surfaces will be increased and spectral contrast lowered 
due to multiple reflections. 

Finally, eq. (2.20) must be adjusted for topo- 
graphic effects in order to have predictive capabil- 
ities for rough terrain. To first order, the emitted 
term is unaffected by terrain. Terrain, however, 
shields longwave irradiance from the atmosphere, 
although this is substituted by longwave adjacent 
terrain irradiance, as the surrounding terrain sur- 
face is generally warmer than the sky. These terrain 
influences can be estimated using a DEM and eq. 
(2.8) to modify the previous equation such that 

L?(A) = TT(A) (e,(A) + (1 - e,(A)) 

x {L^X)iT\\-V f ) + 2itL\{T t ,X)V f ^ 

+ L\{\), ( 2 . 21 ) 

where L\ is longwave irradiance from the adjacent 
terrain and T, is the temperature of the adjacent 
terrain. Here it is assumed that the adjacent terrain 
and sky are isothermal. These assumptions are un- 
realistic, the terms multiplying (1 — e s (A)) are better 
regarded as spatial averages over the sky hemi- 
sphere, and T, is better regarded as a variable, 
especially for work on snowfields and glaciers, 
where the terrain may be bare and have a strongly 
different temperature than the snow. Parameter Vf 


can be called a “form factor” in radiosity theory, 
where it describes the solid angle subtended by one 
surface element viewed from another. 

For many snow or ice studies at high altitude, the 
radiosity effect is relatively small because the reflec- 
tivity of the snow is low in the longwave part of the 
spectrum. For example, if the snow surface is at 273 
K and the adjacent bare hillsides (e s = 0.95) are at 
303 K, and 10% of the sky hemisphere is shielded, 
the apparent temperature of the snow surface is 
increased only ~ 0.2 K. Only if the rocks are warmer 
or substitute more of the sky hemisphere is this a 
measurable effect in imagery (e.g., ASTER). 

It is noteworthy that the atmospheric terms 
and Lj, are controlled by the vertical distributions of 
gases and temperature, with water vapor playing a 
key role. Parameter I.' is heavily weighted by con- 
tributions from the lower atmosphere, but L' p is 
strongly influenced by absorption due to gases 
throughout the atmospheric column, a process that 
is independent of temperature. 

For terrain with low relief, once atmospheric 
variables and emissivities are known, T s is esti- 
mated according to eq. (2.19). In areas of high 
relief, however, especially if hillsides adjacent to 
snowfields or glaciers are bare of snow, it may be 
necessary also to know the temperature and emis- 
sivity of adjacent surfaces in order to solve for T s . 
Here it is best to use a temperature/emissivity 
separation algorithm to estimate both temperatures 
and emissivities from ASTER images (Gillespie et 
al. 1998). Although noise-equivalent change in tem- 
perature (evaluated at 300 K) for ASTER is <0.3 
K, temperatures can typically be recovered only 
within ±1 K, due to compounded uncertainties in 
atmospheric variables, calibration coefficients, and 
assumed emissivity. 


2.3 SURFACE-ENERGY INTERACTIONS 
2.3.1 Snow 

Use of optical, infrared, and microwave sensors 
provides synergy that may allow extraction of 
important snowpack properties to improve our 
understanding of matter- energy interactions in a 
snowpack (e.g., Hall et al. 2005). Snow-covered 
areas can generally be detected most accurately 
using optical data, and augmented by microwave 
data when darkness and clouds prevail. Snow 
albedo can be estimated using data from the reflec- 
tive part of the spectrum. Snow water equivalent, 
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Figure 2.8. Reflectance of snow with different effec- 
tive particle size. Reflectivity data obtained from the 
ASTER Spectral Library (http:/ lspeclib.jpl.nasa.gov) , 
which was collected at the Johns Hopkins University 
IR Spectroscopy Lab. The laboratory spectra here may 
be compared with idealized radiative transfer-based 
computed spectra of pure ice shown in Figure 3.4 of 
this book (Furfaro et al.’s chapter). 

however, is an important hydrologic parameter, 
and can be estimated using microwave wavelengths 
that penetrate the snowpack or emanate from the 
ground. 

One of the first observations made in the study of 
early satellite images of the Earth’s surface was 
coverage of some surfaces by high-albedo snow. 
Surface albedo is the ratio of upwelling reflected 
flux to downwelling incident irradiance. New, 
fine-grained snow cover reflects more than ~80% 
of incident solar radiation, while older, coarse- 
grained and/or melted and refrozen snow tends to 
reflect far less radiation (Fig. 2.8), especially when it 
has been covered by soot/dust or volcanic ash. 
Broadband albedo is reflectance across the reflective 
part of the solar spectrum. Broadband albedo 
decreases when grain size increases as the snow ages 
(Choudhury and Chang 1979), and melting causes 
snow grains to grow and bond into clusters (Dozier 
et al. 1981, Warren 1982, Dozier 1989). 

With the onset of surface melting and associated 
grain size increases, near-infrared reflectance 
decreases dramatically (e.g., Fig. 2.8; Warren and 
Wiscombe 1980, Warren 1982). The near-infrared 
albedo of snow is very sensitive to snow grain size 
while visible albedo is less sensitive to grain size, but 
is affected by snow impurities. For example, for the 
case of just 3 ppm of soot in snow, see the computed 
spectra in Figure 3.6 (Furfaro et al.’s chapter), and 
compare also with their computed Figure 3.4 for 


pure snow and our Figure 2.8 for actual snow. 
Warren and Wiscombe (1980) considered the case 
for snow containing aerosol impurities. Effective 
snow grain radii typically range in size from 
~ 50 pm for new snow, to 1 mm for wet snow con- 
sisting of clusters of ice grains (Warren 1982). Snow 
albedo may decrease by >25% within just a few 
days as grain growth occurs (Nolin and Fiang 
2000). The albedo of a snow cover, if it is thin, is 
also influenced by the albedo of the land, especially 
when vegetation extends above the snow surface. 

Snow albedo can be approximated over large 
areas by deriving the linear relationship between 
the brightest snow-covered areas such as arctic tun- 
dra and the darkest snow-covered forest in a scene. 
Albedos can be assigned 0.80 and 0.18, respectively, 
as did Robinson and Kukla (1985) using Defense 
Meteorological Satellite Program (DMSP) imagery 
(0.4— 1.1pm). Finear interpolation was then per- 
formed to derive albedo. Today, snow albedo can 
be estimated directly from satellite data using the 
Moderate-resolution Imaging Spectroradiometer 
(MODIS) (e.g., Klein and Stroeve 2002, Schaaf et 
al. 2002). 

A near-surface global algorithm has been devel- 
oped to map snow albedo using MODIS data. In 
deriving albedo, atmospherically corrected MODIS 
surface reflectance in individual MODIS bands for 
snow-covered pixels located in nonforested areas 
is adjusted for anisotropic scattering effects using 
a discrete ordinates radiative transfer model 
(DISORT) and snow optical properties (Klein 
and Stroeve 2002). Currently in the algorithm, 
snow-covered forests are considered to be Famber- 
tian reflectors. The adjusted spectral albedos are 
then combined into a broadband albedo measure- 
ment using a narrow-to-broadband conversion 
scheme developed specifically for snow (Nolin 
and Fiang 2000, Klein and Stroeve 2002). Thus 
derived, a near-global snow albedo product is 
available as part of the suite of MODIS snow cover 
products and has been since February 2000 (Hall 
and Riggs 2007). 

In the thermal region, eq. (2.21) is underdeter- 
mined, and a solution for T s requires independent 
information, especially for atmospheric conditions. 
For ASTER images, these data have come from 
NCAR/CIRES reanalysis data from radiosonde 
probes of water, temperature, and pressure profiles. 
These profiles are converted to L and L p using the 
MODTRAN radiative transfer model. For work on 
snow and ice, emissivity is generally measured in the 
laboratory or taken from a spectral library. For 


34 Theoretical foundations of remote sensing for glacier assessment and mapping 


snow and ice, terrain irradiance is less important in 
the thermal infrared than in the visible, near-infra- 
red, and shortwave infrared, because snow and ice 
have reflectivities < 5% between wavelengths of 8 
and 12 pm (Fig. 2.7). Nevertheless, steep, bare rock 
walls above glaciers will have a measurable effect, 
all the more important because the reflectivity of 
snow and ice differ at longer wavelengths, with 
snow a factor of two or more less reflective than 
ice. Estimation of snow and ice temperatures is best 
made at ~ 10.5 pm, because the emissivities are 
about the same and no distinction between them 
is necessary in advance. 

2.3.2 Glaciers 

Glacier surfaces are composed of a variety of 
materials that exhibit unique reflectance patterns. 
High spatial variability in reflectance results from 
variations in the debris cover, ice properties, surface 
moisture content, and snow and firn properties. For 
example, in ablation areas, debris cover lithology, 
depth, and surface roughness significantly affect 
ablation and moisture content, while reflectance 
from ice surfaces is dominated by properties such 
as air bubble and englacial load characteristics. 
Furthermore, centimeter-scale roughness from 
cryoconite holes in ice can produce additional spec- 
tral reflectance variation. In the accumulation area, 
reflectance is controlled by the magnitude of meta- 
morphism of ice crystals, firm density, and moisture 
content (Rees 2006). Collectively, surface matter 
and property variation, multiscale topographic 
effects, and variations in solar geometry result in 
anisotropic reflectance (Greuell and de Wildt 1999). 


Greuell and de Wildt (1999) measured and 
empirically modeled the anisotropic reflection of 
the radiation of melting glacier ice on the 
Morteratschgletcher Glacier in Switzerland. They 
found that anisotropy increases with longer wave- 
lengths (Fig. 2.9), with increasing 9 h and with 
decreasing a. Wavelength dependence can be 
explained due to changes in the absorption coeffi- 
cient of ice with wavelength (Grenfell and Perovich 
1981). A higher absorption coefficient produces 
lower surface albedo, which also dictates the shorter 
effective pathlength of photons through the ice. 
Statistically, a shorter pathlength governs the 
number of scattering events, which leads to better 
preservation of the angular distribution of photons 
in multiple scattering events, which is described by 
the scattering phase function f p . The amount of 
anisotropy increases with decreasing 6 t for snow- 
covered glaciers, as forward-scattered photons 
penetrate more deeply for small zenith angles, 
and photons are more evenly distributed (Warren 
et al. 1998). Lower albedo is caused by surface 
impurities and debris that can sometimes increase 
surface meltwater, which subsequently affects the 
absorption coefficient in a positive feedback. 

The relationships between fundamental radio- 
metric parameters and various properties of the 
glacier surface can be made explicit by considering 
the basic physical laws that describe the conserva- 
tion of energy for photons being transported across 
particulate media. As photons reach the glacier 
surface they enter the ice-snow-firn medium and 
are subject, in principle, to the same absorption- 
scattering mechanisms that occur in the atmo- 
sphere. Nevertheless, glacier surface properties are 
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Figure 2.9. Broadband BRDFs measured over melting glacier ice: (left) Thematic Mapper band 2, 0, = 50°, 
a = 0.52; (right) Thematic Mapper band 4, 9 = 50°, a = 0.50 (modified after Greuell and de Wildt 1 999). 


Surface-energy interactions 35 


very different than atmospheric properties. Using 
first principles, physical models can be derived to 
compute the spectral radiance leaving the glacier 
surface, as well as to compute the BRDF/BRF 
parameters to provide knowledge of surface behav- 
ior. In the specific case of glaciers, the theoretical 
computation of the BRDF/BRF is generally 
accomplished by assuming the surface to be a 
homogeneous, plane-parallel particulate layer of 
ice or snow, and solving the corresponding radia- 
tion transfer equation (RTE). If the optical proper- 
ties (scattering and absorption coefficient and f p ) of 
ice or snow are known, the conservation of photons 
in the appropriate phase space yields the linearized 
Boltzmann equation such that 

d 

<j>) + L(t, n, <j>) 

= ^J J f p {cos &)L{t, <f>) drf d<jJ . (2-22) 

Here we assume that radiance is a function of only 
one spatial dimension (depth) as well as a function 
of inclination (ji = cos #,) and azimuthal angle (</>,). 
Optical thickness (r) is related to thickness ( h ) via 
the extinction coefficient a ext (r = f[! cr ext (h) dh ). 
Single-scattering albedo, to, is the ratio between 
scattering and extinction, while f p describes the 
probability that a photon moving in direction H 0 
is scattered in the direction H. Moreover,^, depends 
on 0 = fl 0 • H (i.e., the angle between incoming and 
scattered photons). 

Eq. (2.22) must be equipped with appropriate 
boundary conditions and solved numerically to 
compute L. Various numerical techniques are 
available to solve the RTE and include the iterative 
solution of the Ambartsumiam’s integral equation 
(Mishchenko et al. 1999), adding-doubling tech- 
nique (Hansen and Travis 1974), converged discrete 
ordinates method (Ganapol and Furfaro, 2008), 
and analytical discrete ordinates method (Siewert 
2000). Areal and intimate mixtures of snow, ice, 
and debris can be modeled using radiative transfer 
theory (see Chapter 3). 

The computed radiance can be directly used to 
evaluate the theoretical BRDF and/or BRF for 
different types of glaciers. Such parameters can be 
employed to study the effect of particle size, shape, 
and orientation of ice crystals on the reflected sig- 
nal. Examples of a RT-computed BRF for ice par- 
ticle models are shown in Fig. 2.10. Mishchenko et 
al. (1999) considered three ice models including: 
(1) particles with a highly irregular, random fractal 


shape (i.e., snow), (2) homogeneous ice spheres, and 
(3) regular hexagonal ice crystals. Reflected inten- 
sities are displayed in polar plots: each column 
corresponds to an ice model (1,2, and 3); the rows 
show BRF as a function of incident angle (the stars 
on the plots indicate the inclination of the incident 
photon beam). Clearly, the scattering phase func- 
tion for the particle type has a marked effect on the 
BRF. While small differences are perceived when 
looking at the nadir direction across models, large 
variations exist in other directions. Note that for 
model 1, the BRF has less structure/features than 
model 3 (i.e., degree of anisotropy). Indeed, a reg- 
ular hexagonal ice model has a more structured 
phase function than a randomly oriented, fractal- 
like model. 

2.3.3 Water 

Increases in temperature and high ablation produce 
significant quantities of meltwater that eventually 
can form supraglacial and proglacial lakes. The 
reader is directed to the following chapter by Fur- 
faro et al. for further theoretical treatment of this 
topic of radiative transfer in glacial waters. The 
interaction between light and water is dictated by 
the same first principles describing photon propaga- 
tion in particulate media. The physics of interaction 
are described by eq. (2.22) which formalizes the 
conservation of photons either absorbed or scat- 
tered or eventually reemitted in any direction. 
Naturally, there are key differences and similarities 
that markedly distinguish the transport of light in 
water from any other medium. For example, in the 
case of water-light interaction, the visible part of 
the spectrum, 400-700 nm, is of interest because it is 
related to biological activity, and is the region 
where optical transmission is the largest. From a 
modeling point of view, glacier meltwater lakes can 
be modeled using the so-called plane-parallel 
approximation. Here it is generally assumed that 
the water column is layered and incident radiation 
is uniform over the surface. 

Optical properties play a critical role in determin- 
ing the behavior of reflected light as collected by 
observing sensors. Such properties (i.e., absorption, 
refraction, and scattering) are generally difficult to 
characterize. Water is the primary compound that 
also potentially incorporates organic and inorganic 
matter. These constituents exhibit high spatio- 
temporal variability in concentration. Multiple 
components are generally present within the water, 
and it is assumed that the absorption and scattering 
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Figure 2.1 0. Two-angle reflected intensity as a function of inclination ( fi = cos 9 V ) and azimuth (<fi v ) for three 
different snow models. The yellow star indicates the incident solar zenith angle (from Mishchenko et al. 1999). 
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coefficients of individual, optically active constitu- 
ents are additive. We can distinguish between clear 
water and turbid meltwater. For water, optical 
properties are primarily determined by the concen- 
tration of suspended detritus. Indeed, there is a 
formal correlation between absorption ( a p ) and 
scattering (b p ) coefficients and concentration of par- 
ticles (C) which can be described as follows (Prieur 
and Sathyendranath 1981): 

a p ( A) = 0.06^(A)C°' 602 , (2.23) 

C 0 - 62 . (2.24) 

A' p and B p range between 0.12 and 0.45 m _I . For 
turbid water, the optical properties (including spec- 
tral variation in reflectance) are strongly influenced 
by mineral particle concentration. Suspended sedi- 
ments/particles tend to be scatterers rather than 
absorbers. Nevertheless, the mineralogical compo- 
sition of the particles can also be a factor especially 
in the case of iron-rich sediments that feature strong 
absorption bands. 

Turbid glacier meltwaters exhibit an extremely 
forward-peaked scattering, with a peak forward 
to backscatter ratio on the order of 10 6 . Scattering 
has been modeled using a forward-peaked phase 
function derived by experimental data (Petzold 
1972) and analytical models such as the 
Fournier-Fourand model (Fournier and Fourand 
1994). 


2.4 COMPLICATIONS 

The previous sections outline the complexity of 
radiation transfer and the relatively large number 
of variables that influence sensor spectral response. 
Information extraction from satellite data is predi- 
cated upon isolating specific variance components 
of the radiation transfer cascade. This frequently 
requires radiometric calibration to account for the 
significant variables that cause spectral variation. In 
practice, accurate information necessary to utilizing 
physical models to reduce unwanted variation 
caused by the atmosphere or the topography is 
not readily available, and Earth scientists typically 
rely upon empirical approaches to reduce or isolate 
spectral variation. 

A classic example for glacier assessment includes 
topographic normalization of satellite imagery. 
Glaciologists rarely account for all the irradiance 
components and viewing geometry related to 


b P ( A) = 0.54(A) 


topography. Instead, empirical procedures such as 
ratioing or principal component analysis is used to 
reduce topographic effects. This is usually effective 
for thematic mapping applications, although quan- 
titative estimation of surface parameters requires 
more detailed insight and treatment of the radiation 
transfer cascade. 


2.5 SPACE-BASED INFORMATION 
EXTRACTION 

Satellite imagery contains a wealth of information 
for assessing the Earth’s cryosphere (Bishop et al. 
2004a). Extracting these signals and pertinent 
glacier information can be challenging. New tech- 
niques and analytical approaches are frequently 
required to identify and transform spectral varia- 
tion into thematic and quantitative parameters. We 
focus here on the type of information that can be 
extracted from satellite imagery. 

2.5.1 Snow cover 

Snow cover mapping can be accomplished at a 
variety of temporal and spatial scales with high 
accuracy. Today, satellite-borne instruments pro- 
vide detailed information on snow cover in the 
visible, infrared, and microwave parts of the elec- 
tromagnetic spectrum on a daily basis at a global 
scale. In 1997, the Interactive Multi-sensor Snow 
and Ice Mapping System (IMS) began to produce 
operational products daily at a spatial resolution of 
about 25 km, utilizing a variety of satellite data 
(Ramsay 1998). Improvements in the spatial reso- 
lution of the IMS product in February 2004 resulted 
in snow maps at 4 km resolution. In June 1999, 
NOAA/NESDIS ceased producing weekly maps, 
replacing them with the daily IMS product (Ram- 
say 1998). 

Since February 24, 2000 the Moderate Resolu- 
tion Imaging Spectroradiometer (MODIS) sensor 
has been providing daily snow maps at a variety 
of different temporal and spatial resolutions (Hall 
and Riggs 2007). MODIS is a 36-channel imaging 
spectroradiometer that was first launched as part of 
NASA’s Earth Observing System (EOS) program 
on December 18, 1999 on the Terra spacecraft. A 
second MODIS was launched on May 4, 2002 on 
the Aqua spacecraft. Snow maps are available from 
MODIS to serve the needs of local, regional, and 
global modelers and are provided at 500 m, 5 km, 
and 25 km resolution, and as daily, weekly, or 
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monthly composites, including fractional snow 
cover and albedo (Hall and Riggs 2007). Snow 
cover maps are also routinely generated using 
higher resolution satellite imagery coupled with a 
variety of pattern recognition methods. This results 
from the good statistical separability between snow 
reflectance and the reflectance of other land-cover 
classes. 

2.5.2 Ice sheets 

Using Advanced Very High Resolution Radiometer 
(AVHRR) data and, more recently, MOD1S data, it 
is possible to achieve daily mapping of the Green- 
land Ice Sheet. With microwave data, most cloud 
cover is not an issue, but there are also many com- 
plexities in terms of what is being observed, and the 
use of data from optical wavelengths provides a 
wealth of information about ice sheet characteris- 
tics including surface temperature and albedo. 

Comiso (2006) has shown, using monthly 
AVHRR-derived, surface temperature maps, that 
much of the Arctic has been undergoing a steady 
warming since 1981, including a 1.19°C per decade 
increase over Greenland. More recently, Hall et al. 
(2006) provided improved temporal and spatial 
detail using 5 km resolution MODIS surface tem- 
perature maps that show the temperature distribu- 
tion in the various facies of the ice sheet. The glacier 
facies (Benson 1996) are delineated based on snow 
and ice characteristics, and subsequently changes of 
the facies boundaries over time can be correlated 
with ice sheet mass balance. For example, using 
surface temperature maps, an increase in the eleva- 
tion of the boundary between the percolation and 
dry-snow facies could be monitored over time to 
reveal associated temporal changes in ice sheet mass 
balance. 

Satellite data (AVHRR) have also been used to 
map changes in albedo over the Greenland Ice 
Sheet during the spring and summer months (Knap 
and Oerlemans 1996, Nolin and Stroeve 1997, 
Stroeve et al. 1997, Comiso 2006). Recently, Stroeve 
et al. (2006) have used MODIS data to analyze the 
accuracy of daily albedo measurements from 
MODIS of the Greenland ice sheet. 

Surface albedo over the Greenland Ice Sheet has 
also been measured using data acquired by the EOS 
Multiangle Imaging Spectroradiometer (MISR) 
instrument. MISR images the surface using nine 
discrete, fixed angle cameras, one nadir viewing 
and four viewing angles in forward and aftward 
directions along the spacecraft track. Compared 


with in situ measurements at five different sites, 
the surface albedo derived from two different meth- 
ods using MISR data showed good agreement 
(Stroeve and Nolin 2002). Although the atmosphere 
is relatively thin over the ice sheet, atmospheric 
attenuation is significant in the visible and near- 
infrared wavelengths. 


2.5.3 Alpine glacier mapping 

Glaciers are known to be directly and indirectly 
sensitive to climate variations, and a top priority 
of assessing climate change is to monitor and 
understand glacier fluctuations at a variety of scales 
(Haeberli et al. 1998, Bishop et al. 2004a, Kargel et 
al. 2005, Haeberli 2004), and to quantify the errors 
(Hall et al. 2003). This requires information on 
changes in glacier distribution. 

The nature and uncertainty of information 
extracted from satellite imagery is strongly related 
to sensor characteristics. Relatively high spatial 
resolution and good image geometric fidelity is 
required to resolve smaller alpine glaciers and 
glacier features. High spatial resolution also permits 
spatial analysis of glacier surfaces, such that ice flow 
and structural characteristics can be assessed and 
mapped (e.g.. Bishop et al. 1998). In addition, scene 
coverage should be large enough to enable regional 
analysis. 

Spectral resolution controls spectral differentia- 
tion of glacier surfaces from other landscape feat- 
ures. Glacier surface reflectance is governed by ice 
structures and the concentration of air bubbles 
(Winther 1993), variations in grain size of snow/firn 
(Kaab 2005b), surface water content (Konig et al. 
2001), and dust and/or supraglacial debris. 

In general, the near and shortwave infrared 
regions of the spectrum are important for snow 
and ice mapping, and many different satellite 
sensors can be used to map alpine glaciers that have 
minimal debris cover. Atmospheric conditions, 
however, such as cloud cover and aerosol influences 
can affect spectral response. The anisotropic nature 
of surface reflectance is also an issue, as multiscale 
topographic effects can significantly influence inter- 
pretation (Bishop et al. 2004a, b). Existing spectral 
bandwidth characteristics, however, make it diffi- 
cult to map debris-covered glaciers as spectral dif- 
ferentiation between rocks and sediment on and off 
glaciers is relatively poor (Bishop et al. 2001, 2004b, 
Paul et al. 2004). Other issues include radiometric 
resolution and sensor gain which can lead to sensor 
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saturation over snow and ice, and repeat coverage 
which affects change detection studies. 

Alpine glacier mapping can be accomplished 
using a variety of methods, each with advantages 
and disadvantages. The most common approach 
used worldwide is based on cursor tracking, which 
is commonly known as on-screen manual digitiza- 
tion from a false-color composite (FCC), and/or 
spectrally enhanced imagery (i.e., ratios or principal 
components). This approach takes advantage of 
spectral variations of landscape features; it is labor 
intensive and fraught with difficulties, as results are 
necessarily dependent upon geographic knowledge 
and experience of the analyst. Glacier boundaries 
can easily be misinterpreted. Other techniques for 
glacier boundary delineation include the image 
ratio approach (e.g., Hall et al. 1989, Paul et al. 
2002), supervised and unsupervised classification 
(e.g., Gratton et al. 1990, Aniya et al. 1996), and 
the use of artificial neural networks (e.g.. Bishop et 
al. 1999). Most of these techniques require image 
preprocessing and work reasonably well for clean 
ice, but they tend to fail for debris-covered glaciers. 
Details regarding mapping glaciers using satellite 
spectral data are covered in Chapter 4. 

Change detection studies have produced esti- 
mates of retreat rates and changing mass balances 
for alpine glaciers around the world. The crucial 
factors for temporal analysis include selecting 
suitable images, identifying change categories, and 
using appropriate change detection algorithms (Lu 
et al. 2004). Space-based glacier change detection 
studies include monitoring of variations in terminus 
position (retreat, advance, and/or surge), down- 
wasting, snow line, debris cover, exposed glacial 
area, supraglacial lake fluctuations, and mass bal- 
ance. Some of these parameters may be interrelated, 
such as glacier terminus position and mass balance 
changes (Haeberli and Hoelzle 1995), and changes 
in supraglacial debris thickness and cover and 
glacier surface downwasting. Paul et al. (2004) indi- 
cated that many glaciers are currently downwasting 
instead of retreating. It is important to note, how- 
ever, that methods to measure annual mass balance 
are still being studied (Khalsa et al. 2004). 

2.5.4 Debris-covered glaciers 

Mapping debris-covered glaciers is notoriously dif- 
ficult (Williams et al. 1991, Bishop et al. 2001). 
Different approaches have been attempted, and 
the use of topographic information in addition to 
multispectral imagery is generally required. Paul et 


al. (2004) developed an automated approach to 
map debris-covered glaciers using slope and curva- 
ture derived from a DEM and neighborhood anal- 
ysis. Bishop et al. (2001) used a scale-dependent 
automated object-oriented approach which classi- 
fies debris-covered glaciers and surface features 
such as ice cliffs and moraines. These examples 
demonstrate the feasibility of automated mapping 
of debris-covered glaciers, although the results are 
highly dependent upon the quality and resolution of 
the DEM, and on the sophistication of the 
approach. Research indicates that integrated data 
fusion approaches can be well suited to effectively 
make use of spectral data and multiple forms of 
topographic information. Details regarding DEM 
generation and geomorphometry are covered in 
Chapter 5. 

Data on the supraglacial debris cover character- 
istics of glaciers is sorely needed, as debris cover 
exerts a tremendous influence on the ablation pro- 
cess (Nakawo and Rana 1999). Previous research 
indicates a complex relationship between debris 
cover thickness and ablation, where debris cover 
can either retard or enhance ablation rates depend- 
ing on the spatial variability in debris lithology and 
thickness (Loomis 1970, Mattson et al. 1993). 
Remote-sensing studies indicate that debris cover 
thickness can be highly variable (Bishop et al. 
1995), and field research indicates a general pattern 
of greater debris depths towards the terminus and 
lateral edges of many alpine glaciers (e.g., Gardner 
and Jones 1993). These patterns are caused by mass 
movement and ice flow dynamics that transport 
debris to the terminus and edges via horizontal 
and vertical advection. Increases in surface ablation 
and glacier surface downwasting overtime, and 
sometimes increases in delivery of debris to the 
glacier by landslides and other mass movements, 
can result in increases in the percentage of debris 
cover. Conversely, acceleration of flow due to 
positive mass balance or surging sometimes can 
move debris to terminal and lateral moraines and 
partially clear away debris from much of the glacier 
surface. Consequently, field or remote-sensing 
assessment of debris cover patterns and debris 
depth variability is essential for estimating 
ablation, glacier mass balance, glacier sensitivity 
to climate change, and water supply from the abla- 
tion area. 

Satellite imagery can be used to map variations in 
supraglacial lithologies. For example, Nakawo et 
al. (1993) were able to differentiate and map areas 
of granitic and schistose debris over the Khumbu 
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Glacier in Nepal. In other areas in the Himalaya, 
the predominant debris lithologies can be mapped 
and include differentiation of: (1) granitic and 
black metasediment debris in the Hindu Kush, 

(2) granodiorite, weathered granodiorite, and black 
metasediments in the Hunza Himalaya, and 

(3) granite, gneiss, and black and red metasediment 
debris in the Karakoram Himalaya. Such capabil- 
itites also exist for debris-covered glaciers in Alaska 
and the Andes Mountains. Difficulties in automated 
mapping and lithologic characterization of debris 
result from spectral variation caused by glacier 
topography, moisture-laden debris, and supra- 
glacial lakes. This is especially the case when debris 
lithology is predominately mafic in composition 
and numerous features exhibit relatively low reflec- 
tance throughout the spectrum. Whereas finely 
intermixed debris and ice (areal or intimate mix- 
tures) can be assessed quantitatively by analysis 
of satellite imagery, and areas of continuous debris 
can be mapped, debris thickness is very difficult to 
assess by remote sensing. 

Research regarding the estimation of ablation 
under a debris layer indicates that the surface tem- 
perature of debris may be used to estimate thermal 
resistance (Nakawa and Young 1981, Nakawo et al. 
1993). This may work accurately only if the debris is 
relatively thin, as it sometimes is. However, some 
efforts to assess debris thicknesses, or to assess the 
edges of glaciers using thermal data over debris- 
covered glaciers, may be thwarted or rendered par- 
tially effective if the debris cover is thick or is of 
uncertain grain size. Grain size and lithology- 
dependent thermal conductivities may vary signifi- 
cantly, causing uncertainty in the results of these 
approaches. More troubling, if debris is more than 
a few centimeters thick, thermal conductive equilib- 
rium is not obtained; rather, temperatures can be 
controlled by the thermal inertia of the debris 
(which includes heat capacity as well as thermal 
conductivity), basically causing internal storage 
and release of heat in the debris. Current models 
inadequately address this problem. 

Thermal inertia variations are primarily caused 
by variations in thermal conductivity which is 
primarily controlled by particle size distribution. 
A surface composed of large rock boulders will 
have a relatively high thermal inertia, while smaller 
particles and fine-grained dust will have lower ther- 
mal inertia. Diurnal temperature variations will not 
be influenced by debris conditions deeper than one 
diurnal thermal skin depth. Debris depth tempera- 
ture profiles have been documented to exhibit non- 


linear variation, presumably due to thermal inertia. 
Moisture conditions, including evaporative cooling 
on debris surfaces, also must be accounted for, but 
there is no easy way to do this using remote-sensing 
data. Though thermal data provide information 
related to glacier mapping, great care must be exer- 
cised when using satellite temperature estimates to 
assess debris depths. 

2.5.5 Snow line and ELA 

The transient snow line (TSL) is the sometimes 
irregular edge between snow and bare glacier ice 
on a glacier. Snow line elevation is strongly con- 
trolled by latitude, topography, and regional 
climatic systems that govern irradiant flux, precipi- 
tation, and snow accumulation. Seasonal variations 
in the altitude of the TSL are significant. It can be 
estimated using satellite imagery and a DEM. The 
highest altitude of the TSL during the ablation 
season is often considered to be equivalent to the 
glacier equilibrium line altitude (ELA), the demar- 
cation between a glacier’s accumulation and abla- 
tion zones (Carrivick and Brewer 2004, Kulkarni et 
al. 2004). 

Given the paucity of regional mass balance infor- 
mation, there is an urgent need to produce estimates 
of mass balance for individual glaciers over exten- 
sive regions. Space-based assessments of the TSL 
permit this, as field-based mass balance estimates 
can be calibrated with the accumulation area ratio 
(AAR) (Kulkarni et al. 2004). The AAR is the ratio 
between the accumulation area and total glacier 
area. When mass balance data exist for several 
glaciers, the functional form of the relationship 
between the AAR and the specific mass balance 
can be used within a region, given numerous 
assumptions (Khalsa et al. 2004). 

For example, Kulkarni et al. (2004) estimated the 
specific mass balance of 19 glaciers in the Himachal 
Pradesh region of the India Himalaya, based upon 
glaciological mass balance data obtained from 1982 
to 1988. A strong regression relationship was 
obtained and used with the LISS-III sensor of the 
Indian Remote Sensing Satellite (Fig. 2.11). It is 
important to recognize that significant uncertainties 
exist and are related to: (1) inability of the glacio- 
logical method to adequately account for the spatial 
variation in ablation and accumulation; (2) vari- 
ability in topographic conditions that govern short- 
wave and longwave irradiance fluxes, and therefore 
ablation and ELA variation; and (3) spatial varia- 
tions in accumulation rates due to variations in the 
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Figure 2.1 1 . Relationship between the accumulation 
area ratio and the specific mass balance for Shaune 
Garang and Gor Garang glaciers in the India Himalaya. 
Reproduced from data in Kulkarni et al. (2004). 


influence and/or coupling of regional climate sys- 
tems. In addition, glacier hypsometry and regional 
topographic variations must also be considered. 

Nonetheless, a variety of statistical methods have 
been developed and tested to estimate ELAs using 
reflectance and topographic data. The most com- 
mon include the AAR approach, the area-altitude 
balance ratio (AABRs), the area-altitude balance 
index (AABI), and median elevation of a glacier 
(MEG) (Meierding 1982, Furbish and Andrews 
1984, Kulkarni 1992, Benn and Evans 1998, Osmas- 
ton 2005). Other popular methods to reconstruct 
previous ELAs include maximum elevation of 
lateral moraines (MELM) and various toe-to-head- 
wall altitude ratios (THAR) (Nesje 1992, Torsnes et 
al. 1993). Several studies carried out for snow line 
assessment from the last glacial maximum (LGM) 
using these methods were reviewed by Mark et al. 
(2005). 

2.5.6 Ice flow velocities 

The velocity field of a glacier is an important 
glaciological parameter that governs a variety of 
processes, and controls glacier geometry and extent. 
Ice velocities can be highly variable, and diurnal 
and annual fluctuations have been observed (e.g., 
Jansson 1996, Nienow et al. 2005). Microwave data 
and optical imagery are commonly used to estimate 
ice flow velocities. With the latter, multitemporal 
satellite imagery permits the generation of an aver- 
age horizontal ice flow velocity field over the period 
of separation of the images. 


A basic technique for extracting the velocity field 
is based on feature tracking, where identifiable 
features on the glacier (e.g., crevasses or snow 
dunes) are tracked in images acquired at different 
times. Feature-tracking research has focused pri- 
marily on large polar glaciers and ice streams (Kaab 
2005a, Scambos et al. 1992), although this tech- 
nique can be applied to smaller alpine glaciers 
(Kaab 2005a, b). It works particularly well for 
glacier surfaces where the ice surface structure 
remains unchanged over time. In general, after pre- 
processing the imagery to enhance moving features, 
a computer program searches for correlation 
matches at many points in the image and, for each 
one, produces a subpixel estimate of the horizontal 
velocity based upon the correlation field. Rapidly 
changing alpine glacier surfaces can potentially 
pose a challenge to this approach in terms of uni- 
form velocity fields, as surface features can change 
dramatically due to rapid mass movement and 
transport of debris, ablation changing debris cover 
conditions and ice structures, supraglacial lake 
development, evolution, and drainage, and non- 
linear spatial variations in ice flow direction. 
Furthermore, spectral saturation and heavy debris 
cover can limit the production of usable correla- 
tions. Consequently, more research into developing 
new methods to account for supraglacial reflectance 
variability is warranted. Details regarding feature- 
tracking methodologies are covered in Chapters 4 
and 5. 

Kaab (2005a) identified the need for developing a 
set of comprehensive goals and strategies for inven- 
torying ice flow dynamics. Ice flow is a fundamental 
process of glaciers, and the spatial and temporal 
variability in ice flow fields can provide us with 
new insights into understanding the relationships 
between climate, glaciers, and mountain geo- 
dynamics. For example, Kaab (2005a) used ASTER 
data to estimate glacier flow velocities in the Bhutan 
Himalaya. Alpine glaciers in the Himalaya may be 
classified as summer or winter accumulation type, 
receiving precipitation predominately from the 
monsoon or the westerlies. They receive precipita- 
tion from both, however, and understanding the 
influence of source partitioning and changes in par- 
titioning over time is an important topic that has 
not been adequately addressed. We might expect 
spatial variations in ice flow velocity to be strongly 
associated with ice flux (i.e., snow accumulation). 
An example of this is shown in Fig. 24.6, where 
higher velocities are associated with the glaciers 
on the north side, where basal conditions and 
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orography lead to velocities of up to 200 m yr~’ . On 
the south side, higher ablation and rainfall produces 
thermokarst features, more supraglacial lakes, and 
lower velocities. 


2.6 NUMERICAL MODELING 

The integration of remote sensing-derived param- 
eters and numerical modeling which has been pre- 
viously discussed represents a more sophisticated 
approach to the study of glacier dynamics and 
the interpretation of causal mechanisms. Informa- 
tion generated from satellite imagery and DEMs 
can be used to establish boundary conditions, esti- 
mate input parameters, and constrain simulations. 
We briefly discuss the significance of numerical 
models for understanding glacier fluctuations. 


2.6.1 Climate modeling 


While mapping and monitoring the present con- 
ditions of ice sheets and alpine glaciers is a funda- 
mental and necessary task, predictions of future 
mass balances must rely on knowledge of what 
the climate will likely be sometime in the future. 
Such knowledge can only come from numerical 
modeling of the climate system. In addition, general 
circulation models (GCMs) can be used for paleo- 
climate studies in which direct comparisons with 
proxy climate data recovered from various regions 
may be made. 

The equations that are used to describe the 
evolution of the atmosphere are based on the 
hydrodynamic equations of motion combined with 
thermodynamic laws. They conserve mass, momen- 
tum, and energy. The momentum equations that 
express the three-dimensional force balance may 
be written as follows in the coordinates of x (longi- 
tude), y (latitude), and p (pressure, which is used as 
a vertical coordinate rather than height z): 
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where U, V, and W are the zonal, meridional, and 
vertical wind speeds, respectively; p a is atmospheric 
density; p is air pressure at latitude <p; Q is the 
rotation rate of the Earth; and F x , F y , and F- are 
frictional terms whose form varies depending on 
what processes are resolved by the model. Sub- 
grid-scale motions (i.e., unresolved motion) will 
be parameterized in some fashion and included in 
these terms. On the left-hand side of eqs. (2.25), 

(2.26) , (2.27), the second, third, and fourth terms 
represent three-dimensional advection. The first 
term on the right-hand side is the acceleration pro- 
duced by pressure gradient forces, followed by 
Coriolis acceleration terms. Gravitational accelera- 
tion, g, is seen in the vertical momentum — eq. 

(2.27) . For synoptic-scale atmospheric flow it is 
quite common to simplify and use the hydrostatic 
approximation, which retains only two of the 
terms: 


| Z=~Pa9- (2-28) 

The hydrostatic approximation forces the vertical 
pressure gradient to be balanced by buoyancy. This 
approximation is good for atmospheric motion on 
the global scale, and means that the full set of 
equations is easier to solve numerically. However, 
it comes at a price: convection in the atmosphere is 
precluded under the hydrostatic assumption 
because vertical acceleration is assumed to be zero. 
(There can still be vertical motion as a result of mass 
continuity and horizontal convergence/divergence, 
however.) The benefits of assuming hydrostatic 
balance are therefore tempered by the need to 
parameterize convection and cloud formation 
processes in some fashion. 

An equation of state is required and the assump- 
tion that the atmosphere is an ideal gas is quite a 
good one for normal ranges of atmospheric tem- 
peratures and pressures: 

P = PaRTa , (2.29) 

where T a is the temperature of the atmosphere and 
R is the dry-air gas constant. 

Invocation of energy conservation for an ideal 
gas with a specific heat at constant pressure of c p 
and a thermal expansion coefficient a te yields the 
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following prognostic thermodynamic equation: 
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where H is an energy term to which latent, sensible, 
and radiative heating contribute. Mass conserva- 
tion is expressed in these coordinates by 
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In addition to these dynamic/thermodynamic equa- 
tions, an expression for the conservation of water is 
required and is usually written in terms of specific 
humidity, q\ 




where E v is the rate of evaporation and P is the rate 
of precipitation. 

Eqs. (2.25) through (2.32) — using either (2.27) or 
(2.28) — represent seven equations for the seven 
unknowns U, V, W, T a , p, p a , and q. These equa- 
tions, in one form or another, are at the core of all 
modern general circulation models. 

The seemingly innocuous heating term on the 
right-hand side of eq. (2.30) contains a wealth of 
physics wherein contributions from remote sensing 
can be made. Radiative transfer through the 
atmosphere is computed for both shortwave and 
longwave radiation and a particular atmospheric 
composition (which can change in its content of 
aerosols and water vapor). The divergence/conver- 
gence of radiative fluxes produces cooling/heating 
terms through H. Latent heating during precipita- 
tion/convection events is computed according to the 
particular convective parameterization used and 
also contributes to H, as does the sensible heating 
associated with parameterized convective motions. 

Surface temperature is a crucial element and is 
computed from a surface energy balance that takes 
into account sensible and latent heat fluxes, net 
shortwave and net longwave radiation. Implicit in 
the calculation is surface albedo. Some general 
circulation models incorporate land surface para- 
meterization schemes in which the land surface can 
respond to climate forcing and hence albedo can 
change according to land cover dynamics. It is com- 
mon to assume the surface to be a graybody such 


that albedo does not depend on wavelength. His- 
torically, the albedo of large ice masses was fixed at 
a particular value and only snow cover could 
change. It is now known that the albedo of glaciers 
has significant seasonal variation due to meltwater 
and surface debris. Crude parameterizations have 
been developed to address this seasonality, but 
remote sensing is required to obtain atmospheric 
and land surface information that better charac- 
terizes the spatiotemporal complexities of these 
important parameters. 

2.6.2 Energy balance modeling 

Together with accumulation, the overall factor that 
governs the mass balance of a glacier is surface 
energy balance. Surface energy budget is composed 
of a variety of processes that determine the magni- 
tude of ablation (ice mass losses due mainly to 
melting and sublimation). It has been demonstrated 
that the spatial and temporal variability of ablation 
is high in mountain regions given altitude ranges, 
topographic conditions, and debris cover variability 
(Bishop et al. 1995, Arnold et al. 2006). The classi- 
cal method based on positive degree day estimation 
of melting provides a crude proxy to part of the 
energy balance formulation given below, but it is 
directly related only to sensible heat. The degree 
day method implicitly assumes that the other terms 
scale in some way with sensible heat, such that there 
is an empirical relation between degree days and 
melting. This is not strictly true, and can be highly 
inaccurate and misleading, and so the use of the 
positive degree day method is decreasing. 

It is generally acknowledged that three- 
dimensional numerical modeling of the energy 
budget is required to characterize the spatio- 
temporal complexity of numerous processes that 
affect ablation, especially on alpine glaciers. Conse- 
quently, input of spatiotemporal data from remote 
sensing and GIS is needed to reduce the uncer- 
tainties associated with key parameters and 
processes such as albedo, irradiance flux, and 
turbulent fluxes. 

A physically based approach for estimating abla- 
tion requires assessment of energy fluxes. The com- 
ponents of an energy flux can be represented as: 

Qn + Qh + Ql + Qg + Qr + Qm = 0 (2.33) 

where Q N is net radiation, Q H is sensible heat flux, 
Q l is latent heat flux ( Q H and Q L are referred to as 
turbulent heat fluxes), Q G is ground heat flux, Q R is 
the sensible heat flux supplied by rain, and Q M is the 
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energy used for melting snow and ice. Positive 
quantities represent energy gains while negative 
magnitudes indicate energy loss at the surface. 

Ablation rate (M a ) estimates generated from 
energy balance modeling are generally deemed reli- 
able (Hock 2005). Ablation is computed as: 


where p w is the density of liquid water and Lf is the 
latent heat of fusion. 

All of the energy flux components are time vari- 
able. They include momentary fluctuations related 
to weather; diurnal, annual, and longer term cycles; 
and long-term trending related to global climate 
change. Ablation and mass balance assessments 
are most usefully and commonly integrated over 
the annual cycle. However, comparison of winter 
and summer mass balances also provides important 
information about ice throughput and the time- 
scales of glacier dynamical responses to perturba- 
tions. 

2.6.2. 1 A Jet radiation 

The net radiation component represents the balance 
between incoming and outgoing energy absorbed 
and emitted by the surface. It is characterized and 
modeled as shortwave (0.15M.0 pm) and longwave 
(4.0-120 pm) radiation fluxes, and is usually the 
dominant energy source. Net radiation is regulated 
by numerous processes and key parameters that 
represent atmospheric, topographic, and surface 
properties. 

The radiation balance can be computed as: 

Q N = E{\-ot)+l}+l} t +L\ (2.35) 

where E is surface irradiance, i) is total longwave 
sky irradiance, and L) is total emitted longwave 
radiation. 

Albedo a can be estimated using satellite 
imagery. It is also possible to measure albedo in 
the field and empirically model the relationship 
between measured albedo and surface radiance. 
This parameter, however, does vary over time, 
and numerous approaches to modeling albedo have 
been developed (Arendt 1999). 

Longwave radiation can be characterized using 
numerous parameterization schemes. In general, 
incoming longwave irradiance (Lr) can be repre- 
sented as: 

L l = (e a aT A a )V f +{e s aTt){\ - V f ), (2.36) 


where a is the Stefan-Boltzmann constant, e a and e s 
are the emissivities of the air and surrounding 
surface. 

Finally, outgoing longwave radiation can be 
calculated from: 

L) = e s aT* + (l -e s )L l . (2.37) 

Satellite thermal imagery can be used to validate 
and constrain energy balance modeling of surface 
temperature. One of the most challenging terms to 
evaluate analytically, but readily assessed by meas- 
urement, is L' 1 because air temperature and emis- 
sivity vary with altitude, and the effective mean 
temperature of the air mass varies with elevation 
angle measured from the zenith. The simplest 
approach is to approximate L) with a single effec- 
tive air temperature and a single effective emissivity, 
and then use field measurements to modify the 
evaluation. However, it will always be true that 
transient weather conditions other than local sur- 
face temperature, such as humidity and lapse rate, 
will alter downwelling sky irradiance. Like all the 
other terms of the energy balance, it is the inte- 
grated annual effect on ablation (hence, dominated 
in most regions by summer conditions) that is the 
most important. 


2. 6. 2. 2 Turbulent heat fluxes 

Sensible and latent heat fluxes vary due to tem- 
perature and moisture gradients and atmospheric 
turbulence. They do not dominate over longer time 
periods, although latent heat flux can be of major 
importance for short-term ablation rates. 

Modeling these components is usually done using 
the bulk aerodynamic expressions that characterize 
turbulent fluxes as: 

Qh = P a c p C H u(9 : - 0 S ) (2.38) 

Ql = PaL v C E u{q- - q s ) (2.39) 

where L v is the latent heat of evaporation, u is mean 
wind speed, 9 Z and 9 S are mean potential tempera- 
tures, q : and q s are mean specific humidities, and 
C H and C r are the exchange coefficients for heat 
and vapor pressure, respectively. Radiation transfer 
and atmospheric modeling are required to estimate 
spatiotemporal variability in these parameters; 
however, field-based measurements from point 
weather stations are typically used. 
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2.6.23 Ground heat 

For temperate glacial ice, surface temperature is 
essentially zero in the daytime with nocturnal freez- 
ing. For cold glacier ice, the surface layer can be an 
energy sink. Surface temperature is regulated by 
heat conduction from englacial heat flux where 
the temperature gradient is a function of depth, 
advection, and thermal conductivity. 

Sources for this energy component include 
geothermal heat flux, friction caused by basal slid- 
ing, and internal ice deformation. Flux across the 
surface can be represented as 

r z ■ dT 

Qg = J PiC l p-^-dz, (2.40) 

where dTjdt is the rate of change in ice tempera- 
ture, pi is ice density, and c‘ p is the specific heat 
capacity of ice. 

Geothermal flux is not large, except in excep- 
tional circumstances (such as glaciers on active vol- 
canoes), as can be appreciated by the magnitudes of 
global mean geothermal heat flux of about 0.08 W 
m 2 . If there exists a net mass balance driven 
mainly by a net loss or gain of 1% of solar constant 
energy flux (i.e., 13 W m~ 2 ), this is over two orders 
of magnitude greater than geothermal flux, which 
therefore can be ignored except in cases of (1) poly- 
thermal ice sheets in severe polar climates, where 
basal melting due to geothermal flux can be a sig- 
nificant fraction of total mass balance, and (2) the 
aforementioned exceptional cases of glaciers on vol- 
canoes, where geothermal flux can be up to three 
orders of magnitude greater than the global mean. 

A further way to see the relative insignificance of 
geothermal flux is to consider that each year 
geothermal global mean heat flux can melt only 
about 8 mm of ice, which is of the order of 1 % of 
the typical negative balances that exist for many 
glaciers. The frictional energy dissipation due to 
viscous deformation of ice and basal sliding is also 
generally a small term, but it can verge toward 
significance, as can be assessed from the total 
change in gravitational potential energy of a col- 
umn of glacier ice as it moves downslope. Since 
glaciers move downslope without rapid accelera- 
tion, it must be that nearly all the gravitational 
potential energy is converted into frictional energy 
and damage to the ice crystal lattice, which ulti- 
mately is also converted to heat as the ice anneals 
its damage. For instance, a vertical parcel of a 
glacier 300 m thick moving downslope and losing 
10 m elevation per year will lose potential energy 


per unit area at a rate of (300 m * 900 kg m~ 3 ) * (9.8 
m s~ 2 ) * (10 m)/(3.16 * 10 7 s/yr) equal to 0.83 W/m 2 
(i.e., about 10 times global mean geothermal flux). 
If the energy is completely dissipated at the bed by 
sliding (as an approximation to glaciers exhibiting a 
high degree of plug flow), this frictional energy is 
sufficient to melt 79 kgm~ 2 of ice per year, or about 
9 cm of ice per year. This ice is melted at the bed, 
but presumably either infiltrates beneath the bed or 
is eventually expelled by basal fluid flow. This is 
characteristically 10% or so of the negative bal- 
ances recorded for many glaciers around the world. 
So it is perhaps not large, but it can potentially be a 
term that is not completely insignificant, yet it is 
normally not dealt with. Changes in this energy 
balance term, however, can pose a stabilizing nega- 
tive feedback because, as glaciers thicken and flow 
more rapidly, they dissipate more energy and melt 
more ice; whereas, when they thin and flow more 
slowly, they dissipate less energy and hence melt less 
ice. 


2.6.3 Glacier mass balance modeling 

Estimating glacier mass balance is complicated 
because ice mass gain or loss can take place at 
the glacier surface, within the glacier, or at the base 
of the glacier. At the base, loss occurs as a result of 
geothermal heat and sliding friction, although sur- 
face loss is most significant. Internally, mass can 
accumulate as meltwater refreezes and losses occur 
due to thermal erosion. Surface accumulation 
includes snowfall, wind drift, avalanches, deposi- 
tion from vapor, and condensation. Surface abla- 
tion includes melting, sublimation, ice avalanches, 
and snow drift. 

Mass balance estimates have been traditionally 
generated using the glaciological method, which 
involves inferring the mass balance from point 
measurements taken from ablation stakes and snow 
pits. Furthermore, energy balance modeling has 
demonstrated the significant variability in ablation 
caused by glacier topography and localized topo- 
graphic parameters that govern shortwave and 
longwave energy fluxes. For alpine glaciers, topo- 
graphic variation is considerable, and relief and 
topographic shielding significantly vary over rela- 
tively short distances. Consequently, the spatial 
complexity in ablation highlights the difficulty of 
using systematic stake measurements for accurately 
inferring mass balance gradients (Arnold et al. 
2006). 
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Remote sensing and GIS permit the use of the 
geodetic method, wherein changes in glacier topog- 
raphy over a given period of time can be used to 
estimate volume change. Satellite-derived DEMs 
can be used to accomplish this. The total mass bal- 
ance of a glacier (B) is: 

B = Ji J b(x,y) dA s dt = p,A V,. (2.41) 

where A s is the surface area of the glacier, b is the 
mass balance at a particular point (x,y), and dV t is 
the ice volume change. This approach can be used 
in a GIS to compute the mean specific mass balance 
(b) which is useful for comparing glaciers. It is 
computed as: 

£ 'Ex 1 EyMx,y)AV,(x,y) 
b = -j ■ (2.42) 

Berthier et al. (2007) used this approach to estimate 
glacier mass balances in the western Himalaya of 
India. They found that remote-sensing estimates 
generally agreed with mass balances estimated in 
the field. The geodetic method, however, has its 
limitations including DEM accuracy and effective 
removal of biases (Nuth and Kaab 2011, Kaab et al. 
2012 ). 

Volume estimates also need to be corrected for 
the density of surface materials. Generation proce- 
dures for DEMs permit altitude accuracy in moun- 
tain environments on the order of 15-30m (Eckert 
et al. 2005, Fujisada et al. 2005). Density correc- 
tions below the equilibrium line altitude (ELA) are 
reasonably accurate (^830-917 kg m -3 ), although 
above the ELA density can vary from 50kgm -3 
to 830 kgm -3 . Finally, the point mass balance, b, 
cannot be determined from the geodetic method 
because glacier elevation is also a function of mass 
flux divergence and convergence due to ice flow. 
Consequently, numerical flow modeling must also 
be used to relate results from the two methods in 
detail. 

Hubbard et al. (2000) described the use of the 
mass continuity equation to compute mass balance 
such that: 

b(x,y) = dh '^ } ^ + V ■ [H,(x,y) v(x,y)\, (2.43) 

where /?,■ is ice elevation of the glacier surface, H, is 
ice thickness, v is the vertically averaged ice velocity 
vector, and V is the horizontal divergence operator. 
The term, dh jdt, can be determined from the geo- 
detic method, while the local mass flux term can be 
derived from numerical three-dimensional ice flow 


modeling or surface ice velocity fields generated 
using feature tracking. Their results demonstrated 
that this approach has potential for deriving high 
spatial and temporal resolution estimates of mass 
balance, the main limitation being in the accurate 
characterization of basal sliding. 

2.7 CONCLUSIONS 

Remote-sensing science and technology are indis- 
pensable for assessing ice sheets and alpine glaciers, 
as thematic and quantitative estimates of glacier 
parameters can be derived from satellite imagery 
in the reflective and thermal regions of the spec- 
trum. Space-based, change detection studies based 
upon spectral and topographic information are 
providing estimates of glacier fluctuations around 
the world. Unfortunately, systematic inventories of 
other important glacier parameters, such as supra- 
glacial debris characteristics, lakes, glacier features, 
and velocity fields have yet to be effectively realized 
due to a variety of factors. Image information 
content, the complexity of the radiation transfer 
cascade, and new analytical approaches provide 
the basis for new opportunities for information 
extraction. These include improved mapping 
capabilities using reflectance anisotropy, scale- 
dependent spatial analysis, and data fusion for 
characterization of glacier surfaces. Existing work, 
coupled with new developments, will provide us 
with more empirical data that will provide new 
insights into explaining glacier fluctuations. 

Interpreting the causal mechanisms responsible 
for glacier fluctuations requires more sophisticated 
analysis that integrates remote-sensing science and 
technology with deterministic modeling. Climate, 
surface energy budget, ice flow, and glacier model- 
ing are required in order to assess glacier dynamics 
and mass balance. Only in this way can a multitude 
of internal and external forcing factors and 
interrelated processes be evaluated objectively to 
quantitatively ascertain glacier response to climate 
change. We can then generate more reliable esti- 
mates of regional mass balance trends and improve 
estimates of the contribution of melting ice sheets 
and alpine glaciers to rising sea level. 
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2.9 NOTATION 

a p Water absorption coefficient [dimensionless] 
A e Area of the solar disk [km 2 ] 

A s Surface area of a glacier [km 2 ] 

A p Pixel area [m 2 ] 

b Mass balance of a glacier at a particular point 
(x,y), water equivalent [mm] 

b Mean specific mass balance, water equivalent 
[mm] 

b p Water scattering coefficient [dimensionless] 

B Total mass balance of a glacier, water 
equivalent [mm] 

c p Atmospheric specific heat at constant 

pressure [J kg -1 K -1 ] 

c'p Specific heat capacity of ice [J kg -1 K 1 ] 

C Concentration of particles in water [mg L -1 ] 
C E Vapor-pressure exchange coefficient 

C H Heat exchange coefficient 

C\ First radiation constant, 3.74151 x 10 8 

[W m -2 pm 4 ] 

C 2 Second radiation constant, 1.43879 x 10 4 

[pm K] 

d Distance [m] 

d es Earth-Sun distance [km] 

E° Exoatmospheric spectral irradiance 

[W m -2 pm -1 ] 

E Surface spectral irradiance [W m -2 pm -1 ] 

E a Diffuse skylight aerosol-scattered spectral 

irradiance [W rrf 2 pm -1 ] 

E, Adjacent terrain spectral irradiance 

[W m -2 pm -1 ] 

E h Direct beam spectral irradiance 

[W m -2 pm -1 ] 

E d Diffuse skylight spectral irradiance 

[W m -2 pm -1 ] 

E g Diffuse skylight ground-backscattered 

spectral irradiance [W m -2 pm -1 ] 


E,. Diffuse skylight Rayleigh-scattered spectral 
irradiance [W m -2 pm -1 ] 

E, Total surface irradiance [W m 2 ] 

E„ Evaporation rate [kg day -1 ] 
f p Scattering phase function [dimensionless] 

f r Bidirectional reflectance distribution function 

[sr ‘] 

F x Atmospheric frictional term for the x 
direction 

F y Atmospheric frictional term for the y 
direction 

F- Atmospheric frictional term for the r 
direction 

g Gravitational acceleration [m s -2 ] 

h Thickness [m] 

hj Glacier ice altitude [m] 

H Atmospheric heating term [W m 2 ] 

77, Ice thickness [m] 

H d Debris depth [m] 

i Incidence angle of illumination 

k Thermal conductivity [W m -1 K -1 ] 

L Shortwave reflected surface radiance 

[W m -2 pm -1 sr -1 ] 

Ly Longwave surface radiance 

[W m -2 pm -1 sr -1 ] 

L° At-satellite measured shortwave spectral 

radiance [W m 2 pm 1 sr -1 ] 

L? At-satellite measured longwave spectral 

radiance [W m -2 pm -1 sr -1 ] 

Total longwave sky irradiance [W m -2 ] 

L\ Longwave spectral irradiance from adjacent 

terrain [W m -2 pm -1 ] 

l) Total emitted longwave radiation [W m -2 ] 

l) p Upwelling longwave atmospheric path 

radiance [W m -2 pm -1 sr -1 ] 

Lf Latent heat of fusion for ice [J kg -1 ] 

L p Additive path radiance [W m -2 pm -1 sr -1 ] 

L v Latent heat of evaporation [J kg - 1 ] 

M Spectral exitance [W m -2 pm -1 ] 

M a Ablation rate [mm day -1 ] 

p Atmospheric pressure [mb] 

P Precipitation rate [cm yr -1 ] 

cj Specific humidity [kg kg -1 ] 

q s Mean specific humidity at surface [kg kg -1 ] 
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q z Mean specific humidity at altitude z [kg kg 

Q c Conductive heat flux through the debris 
[W m~ 2 ] 

Q n Net radiation flux [W m~ 2 ] 

Q h Sensible heat flux [W m~ 2 ] 

Q l Latent heat flux [W m~ 2 ] 

Q g Ground heat flux [W m 2 ] 

Q r Sensible heat flux from rain [W nT 2 ] 

Q m Melting energy [W m 2 ] 

R Dry-air gas constant [m 3 mb kg~' K 1 1 

R,. Bidirectional reflectance factor 
[dimensionless] 

R, Thermal resistance of debris [nr K W _1 ] 

S Cast shadow binary coefficient 
[dimensionless] 

t Time 

T r . Temperature of the Sun [K] 

T a Temperature of atmosphere [K] 

Tj Temperature of ice [K] 

T s Temperature of surface [K] 

T, Temperature of the adjacent terrain [K] 

T Transmission coefficient [dimensionless] 

T^ Total downward atmospheric transmittance 
coefficient [dimensionless] 

T 1 Total upward atmospheric transmittance 
coefficient [dimensionless] 

T,. Rayleigh atmospheric transmittance 
coefficient [dimensionless] 

T„ Aerosol atmospheric transmittance 
coefficient [dimensionless] 

T 0l Ozone atmospheric transmittance coefficient 
[dimensionless]. 

T ?as Miscellaneous gases atmospheric 

transmittance coefficient [dimensionless] 

Th,o Water vapor atmospheric transmittance 
coefficient [dimensionless] 

u Mean wind speed [m s _l ] 

U Zonal wind speed [m s _1 ] 
v Vertically averaged ice velocity vector 
V Meridional wind speed [m s _l ] 

Vf Hemispherical sky view factor coefficient 
[dimensionless] 

Vj Glacier ice volume [m~ 3 ] 

W Vertical wind speed [m s _1 ] 


x Horizontal cartesian coordinate (longitude 
direction) 

y Horizontal cartesian coordinate (latitude 
direction) 

z Vertical cartesian coordinate (altitude 
direction) 

a Surface albedo [dimensionless] 
a a Surface absorption coefficient [dimensionless] 
a psa Planar surface albedo [dimensionless] 
a ssa Spherical surface albedo [dimensionless] 
a te Atmospheric thermal expansion coefficient 
£j Surface emissivity [dimensionless] 
e a Atmospheric emissivity [dimensionless] 
dj Incident solar zenith angle 

6 max Maximum local horizon angle 
6 r Reflected zenith angle 

6, Terrain slope angle 

6 S Mean potential temperature at surface [K] 

9 Z Mean potential temperature at altitude z [K] 
A Wavelength [pm] 

7T Pi constant 

p Reflectance [dimensionless] 

p a Density of atmosphere [kg m -3 ] 

Pi Density of ice [kg m” 3 ] 

p w Density of water [kg nT ’] 

a Stefan-Boltzmann constant [W m -2 K -4 ] 

t Extinction coefficient 

tp Latitude [radians] 

4>i Incident solar azimuth angle 

(j) r Reflected azimuth angle 

(j>, Slope-aspect angle of the terrain 

H Angular velocity of the Earth [rad s -1 ] 
oj Single-scattering albedo [dimensionless] 

V Horizontal divergence operator 

2.10 REFERENCES 

Aniya, M., Sato, H., Naruse, R., Skvarca, P., Casassa, G. 
(1996) The use of satellite and airborne imagery to 
inventory outlet glaciers of the southern Patagonia 
icefield, South America. Photogrammetric Engineering 
and Remote Sensing , 62(12), 1361-1369. 


References 49 


Arendt, A. (1999) Approaches to modelling the surface 
albedo of a high Arctic glacier. Geografiska Annaler, 
81A(4), 477^188. 

Arnold, N.S., Rees, W.G., Hodson, A.J., Kohler, J. 
(2006) Topographic controls on the surface energy 
balance of a high Arctic valley glacier. Journal of Geo- 
physical Research , 111 , F02011, doi: 10.1029/ 

2005JF000426. 

Bamber, L.L., and Kwok, R. (2004) Remote-sensing tech- 
niques. In: J.L. Bamber and A.J. Payne (Eds.), Mass 
Balance of the Cryosphere: Observations and Modelling 
of Contemporary and Future Changes. Cambridge Uni- 
versity Press, Cambridge, U.K., pp. 59-113. 

Bamber, J.L., and Payne, A.J. (2004) Mass Balance of the 
Cryosphere: Observations and Modelling of Contempor- 
ary and Future Changes. Cambridge University Press, 
Cambridge, U.K. 

Barry, R.G. (2006) The status of research on glaciers and 
global glacier recession: A review. Progress in Physical 
Geography, 30(3), 285-306. 

Benn, D.I., and Evans, D.J.A. (1998) Glaciers and Glacia- 
tion, Arnold, London. 

Benson, C.S. (1996) Stratigraphic Studies in the Snow and 
Firn of the Greenland Ice Sheet (Tech. Rep. 70), U.S. 
Army Cold Regions Research and Engineering 
Laboratory, Flanover, NH. 

Berthier, E., Arnaud, Y., Rajesh, K., Sarfaraz, A., Wag- 
non, P., and Chevallier, P. (2007) Remote sensing esti- 
mates of glacier mass balances in the Himachal 
Pradesh (Western Himalaya, India). Remote Sensing 
of Environment, 108(3), 327-338. 

Bishop, M.P., Shroder, J.F., Jr., and Ward, J.L. (1995) 
SPOT multispectral analysis for producing supra- 
glacial debris-load estimates for Batura Glacier, 
Pakistan. Geocar to International, 10(4), 81-90. 

Bishop, M.P., Shroder, Jr., J.F., Hickman, B.L., and 
Copland, L. (1998) Scale dependent analysis of satellite 
imagery for characterization of glacier surfaces in the 
Karakoram Himalaya. Geomorphology, 21, 217-232. 

Bishop, M.P., Shroder, Jr., J.F., and Hickman, B.L. 
(1999) High resolution satellite imagery and neural 
networks for information extraction in a complex 
mountain environment. Geocarto International, 14(2), 
17-26. 

Bishop, M.P., Bonk, R., Kamp, U., Jr., and Shroder, 
J.F., Jr. (2001) Terrain analysis and data modeling 
for alpine glacier mapping. Polar Geography, 24(4), 
257-276. 

Bishop, M.P., Shroder, J.F., Jr., Bonk, R., and Olsen- 
holler, J. (2002) Geomorphic change in high moun- 
tains: A western Himalayan perspective. Global and 
Planetary Change, 32, 311-329. 

Bishop, M.P., Barry, R.G., Bush, A.B.G., Copland, L., 
Dwyer, J.L., Fountain, A.G., Haeberli, W., Hall, 
D.K., Kaab, A., Kargel, J.S. et al. (2004a). Global 
land-ice measurements from space (GLIMS): Remote 


sensing and gis investigations of the Earth’s cryo- 
sphere. Geocarto International, 19(2), 57-84. 

Bishop, M.P., Colby, J.D., Luvall, J.C., Quattrochi, D., 
and Rickman, D.L. (2004b). Remote-sensing science 
and technology for studying mountain environments. 
In: M.P. Bishop and J.F. Shroder, Jr. (Eds.), Geo- 
graphic Information Science and Mountain Geomor- 
phology, Springer/Praxis, Heidelberg, Germany/ 
Chichester, U.K., pp. 147-187. 

Bush, A.B.G. (2000) A positive feedback mechanism for 
Himalayan glaciation. Quaternary International, 65(6), 
3-13. 

Bush, A.B.G. (2001) Pacific sea surface temperature 
forcing dominates orbital forcing of the early Holocene 
monsoon. Quaternary Research, 55, 25-32. 

Carrivick, J.L., and Brewer, T.R. (2004) Improving local 
estimations and regional trends of glacier Equilibrium 
Line Altitudes. Geografiska Annaler, 86A(1), 67-79. 

Chavez, P.S., Jr. (1996) Image-based atmospheric cor- 
rections: Revisited and improved. Photo gramme trie 
Engineering and Remote Sensing, 62(9), 1025-1036. 

Choudhury, B.J., and Chang, A.T.C. (1979) Two-stream 
theory of reflectance of snow. IEEE Transactions on 
Geoscience and Remote Sensing, GE-17(3), 63-68. 

Comiso, J.C. (2006) Arctic warming signals from satellite 
observations. Weather, 61(3), 70-76. 

Diner, D.J., Asner, G.P., Davies, R., Knyazikhin, Y., 
Muller, J.P., Nolin, A.W., Pinty, B„ Schaaf, C.B., 
and Stroeve, J. (1999) New directions in Earth observ- 
ing: Scientific application of multi-angle remote 
sensing. Bulletin of the American Meteorological 
Society, 80, 2209-2228. 

Dozier, J. (1989) Spectral signature of alpine snow cover 
from the Landsat Thematic Mapper. Remote Sensing 
of the Environment, 28, 9-22. 

Dozier, J., Bruno, J., and Downey, P. (1981) A faster 
solution to the horizon problem. Computers and Geo- 
sciences, 7(1), 145-151. 

Dyurgerov, M.B., and Meier, M.F. (2004) Glaciers and 
the study of climate and sea-level change. In: J.L. 
Bamber and A.J. Payne (Eds.), Mass Balance of the 
Cryosphere: Observations and Modelling of Contempor- 
ary and Future Changes. Cambridge University Press, 
Cambridge, U.K., pp. 579-621. 

Eckert, S., Kellenberger, T., and Itten, K. (2005) 
Accuracy assessment of automatically derived digital 
elevation models from Aster data in mountainous ter- 
rain. International Journal of Remote Sensing, 26(9), 
1943-1957. 

Feng, F.X., Schott, J.R., and Gallagher, T.W. (1993) 
Comparison of methods for generation of absolute 
reflectance factor values for BRDF studies. Applied 
Optics, 32(7), 1234-1242. 

Fournier, G.R., and Fourand, J.L. (1994) Analytic phase 
function for ocean water. In: J.S. Jafle (Ed.), Ocean 
Optics XII. International Society for Optics and Photo- 
nics, Bellingham, WA, pp. 194-201. 


50 Theoretical foundations of remote sensing for glacier assessment and mapping 


Fujisada, H., Bailey, G.B., Kelly, G.G., Hara, S., and 
Abrams, M.J. (2005) ASTER DEM performance. 
IEEE Transactions on Geoscience and Remote Sensing, 
43(12), 2707-2714. 

Furbish, D.J., and Andrews, J.T. (1984) The use of 
hypsometry to indicate long-term stability and 
response of valley glaciers to changes in mass transfer. 
Journal of Glaciology, 30(105), 199-211. 

Ganapol, B.D., and Furfaro, R. (2008) The Art of Analy- 
tical Benchmarking (Lecture Notes in Computational 
Science and Engineering No. 62), Springer, New York, 
pp. 105-134. 

Gardner, J.S., and Jones, N.K. (1993) Sediment transport 
and yield at the Raikot Glacier, Nanga Parbat, Punjab 
Himalaya. In: J.F. Shroder Jr. (Ed.), Himalaya to the 
Sea: Geology, Geomorphology and the Quaternary, 
Routledge, London, pp. 184-197. 

Giles, P.T. (2001) Remote sensing and cast shadows in 
mountainous terrain. Photogrammetric Engineering 
and Remote Sensing, 67(7), 833-839. 

Gillespie, A.R., Cothern, J.S., Matsunaga, T., Roku- 
gawa, S., and Hook, S.J. (1998) Temperature and 
emissivity separation from Advanced Spaceborne 
Thermal Emission and Reflection Radiometer 
(ASTER) images. IEEE Transactions on Geoscience 
and Remote Sensing, 36, 1113-1126. 

Gratton, D.J., Howarth, P.J., and Marceau, D.J. (1990) 
Combining DEM parameters with Landsat MSS and 
TM imagery in a GIS for mountain glacier character- 
ization. IEEE Transactions on Geoscience and Remote 
Sensing, GE-28, 766-769. 

Grenfell, T.C., and Perovich, D.K. (1981) Radiation 
absorption coefficients of polycrystalline ice from 
400-1400 nm. Journal of Geophysical Research, 
86(C8), 7447-7450. 

Greuell, W., and de Wildt, M. (1999) Anisotropic reflec- 
tion by melting glacier ice: Measurements and para- 
metrizations in Landsat TM bands 2 and 4. Remote 
Sensing of Environment, 70, 265-277 . 

Gueymard, C. (1995) SMARTS2, A Simple Model of the 
Atmospheric Radiative Transfer of Sunshine: Algo- 
rithms and Performance Assessment (Tech. Rep. 
FSEC-PF-270-95), Florida Solar Energy Center, 
Cocoa, FL. 

Haeberli, W. (1998) Historical evolution and operational 
aspects of worldwide glacier monitoring. In: W. Hae- 
berli, M. Hoelzle, and S. Suter (Eds.) (1998) Into the 
Second Century of Worldwide Glacier Monitoring: 
Prospects and Strategies, UNESCO, Paris, pp. 35 — 47. 

Haeberli, W. (2004) Glacier and ice caps: Historical back- 
ground and strategies of world-wide monitoring. In: 
J.L. Bamber and A.J. Payne (Eds.), Mass Balance of 
the Cryosphere: Observations and Modelling of Contem- 
porary and Future Changes. Cambridge University 
Press, Cambridge, U.K., pp. 559-578. 


Haeberli, W., and Beniston, M. (1998) Climate change 
and its impacts on glaciers and permafrost in the Alps. 
Ambio, 27(4), 258-265. 

Haeberli, W., and Hoelzle, M. (1995) Application of 
inventory data for estimating characteristics of and 
regional climate-change effects on mountain glaciers: 
A pilot study with the European Alps. Annals of 
Glaciology, 21, 206-212. 

Haeberli, W., Hoelzle, M., and Suter, S. (Eds.) (1998) Into 
the Second Century of Worldwide Glacier Monitoring: 
Prospects and Strategies, UNESCO, Paris. 

Hall, D.K., and Riggs, G.A. (2007) Accuracy assessment 
of the MODIS snow-cover products. Hydrological 
Processes, 21(12), 1534-1547, doi: 10.1002/hyp.6715. 

Hall, D.K., Chang, A.T.C., Foster, J.L., Benson, C.S., 
and Kovalick, W.M. (1989) Comparison of in situ and 
Landsat derived reflectance of Alaskan glaciers. 
Remote Sensing of the Environment , 28, 23-31. 

Hall, D.K., Bayr, K.J., Schoner, W., Bindschadler, R.A., 
and Chien, J.Y.L. (2003) Consideration of the errors 
inherent in mapping historical glacier positions in Aus- 
tria from the ground and space (1893-2001). Remote 
Sensing of Environment, 86(4), 566-577. 

Hall, D.K., Kelly, R.E.J., Foster, J.L., and Chang, 
A.T.C. (2005) Encyclopedia of Hydrologic Sciences, 
John Wiley & Sons, Chichester, U.K. 

Hall, D.K., Williams, R.S., Jr., Casey, K.A., DiGiro- 
lamo, N.E., and Wan, Z. (2006) Satellite-derived, 
melt-season surface temperature of the Greenland Ice 
Sheet (2000-2005) and its relationship to mass balance. 
Geophysical Research Letters, 33, LI 1501, doi: 
10. 1029/2006GL026444. 

Hansen, J.E., and Travis, L.D. (1974) Light scattering 
in planetary atmospheres. Space Science Review, 16, 
527-610. 

Hock, R. (2005) Glacier melt: A review of processes and 
their modelling. Progress in Physical Geography, 29(3), 
362-391. 

Hubbard, A., Willis, I., Sharp, M., Mair, D., Nienerow, 
P., Hubbard, B., and Blatter, H. (2000) Glacier mass- 
balance determination by remote sensing and high- 
resolution modelling. Journal of Glaciology, 46(154), 
491^198. 

Jansson, P. (1996) Dynamics and hydrology of small 
polythermal valley glacier. Geografiska Annaler, 
78A(2/3), 171-180. 

Kaab, A. (2005a) Combination of SRTM3 and repeat 
ASTER data for deriving alpine glacier flow velocities 
in the Bhutan Himalaya. Remote Sensing of Environ- 
ment, 94, 463 — 474. 

Kaab, A. (2005b) Remote Sensing of Mountain Glaciers 
and Permafrost Creep (Vol. 48 of Physical Geography 
Series), Geographisches Institut der Universitat, 
Zurich. 

Kaab, A., Berthier, E., Nuth, C., Gardelle, J., and 
Arnaud, Y. (2012) Contrasting patterns of early 


References 51 


twenty-first-century glacier mass change in the Hima- 
layas. Nature , 488, 495^198. 

Kargel, J.S., Abrams, M.J., Bishop, M.P., Bush, A., 
Hamilton, G., Jiskoot, H., Kaab, A., Kieffer, H.H., 
Lee, E.M., Paul, F. et al. (2005) Multispectral imaging 
contributions to global land ice measurements from 
space. Remote Sensing of the Environment , 99, 187-219. 

Khalsa, S.J.S., Dyurgerov, M.B., Khromova, T., Raup, 
B.H., and Barry, R.G. (2004) Space-based mapping of 
glacier changes using ASTER and GIS tools. IEEE 
Transactions on Geoscience and Remote Sensing, 42, 
2177-2182. 

Klein, A.G., and Stroeve, J. (2002) Development and 
validation of snow albedo algorithm for the MODIS 
instrument. Annals of Glaciology, 34, 45-52. 

Knap, W.H., and Oerlemans, J. (1996) The surface 
albedo of the Greenland Ice Sheet: Satellite-derived 
and in situ measurements in the Sondre Stromfjord 
area during the 1991 melt season. Journal of Glaciol- 
ogy, 42(141), 364-374. 

Konig, M., Winther, J.G., and Isaksson, E. (2001) 
Measuring snow and glacier ice properties from satel- 
lite. Reviews of Geophysics, 39, 1-27. 

Kotlyakov, V.M., Serebrjanny, L.R., and Solomina, 
O.N. (1991) Climate change and glacier fluctuations 
during the last 1,000 years in the southern mountains 
of the USSR. Mountain Research and Development, II, 
1 - 12 . 

Kulkarni, A.V. (1992). Mass balance of Himalayan 
glaciers using AAR and ELA methods. Journal of 
Glaciology, 38, 101-104. 

Kulkarni, A.V., Rathore, B.P., and Alex, S. (2004) 
Monitoring of glacier mass balance in the Baspa basin 
using accumulation area ratio method. Current 
Science, 86(1), 185-190. 

Loomis, S.R. (1970) Morphology and ablation processes 
on glacier ice. Proceedings of the Association of Amer- 
ican Geographers, 2, 88-92. 

Lu, D., Mausel, P., Brondzio, E., and Moran, E. (2004) 
Change detection techniques. International Journal of 
Remote Sensing, 25(12), 2365-2401. 

MacClune, K.L., Fountain, A.G., Kargel, J.S., and 
MacAyeal, D.R. (2003) Glaciers of the McMurdo 
dry valleys: Terrestrial analog for Martian polar sub- 
limation. Journal of Geophysical Research, 108(E4), 
5031-5043. 

Maisch, M. (2000) The longterm signal of climate change 
in the Swiss Alps: Glacier retreat since the end of the 
Little Ice Age and future ice decay scenarios. Geografia 
Fisica e Dinamica Quaternaria, 23(2), 139-151. 

Mark, B.G., Harrison, S.P., Spessa, A., New, M., Evans, 
D.J.A., Helmens, K.F. (2005) Tropical snowline 
changes at the Last Glacial Maximum: A global assess- 
ment. Arctic and Alpine Research , 138/139, 168-201. 

Mattson, L.E., Gardner, J.S., and Young, G.J. (1993) 
Ablation on debris covered glaciers: An example from 
the Rakhiot Glacier, Punjab, Himalaya. In: Snow and 


Ice Hydrology (No. 218), International Association of 
Hydrological Sciences, Rennes, France. 

McClung, D.M., and Armstrong, R.L. (1993) Temperate 
glacier time response from field data. Journal of 
Glaciology, 39(132), 323-326. 

Meier, M.F., and Wahr, J.M. (2002) Sea level is rising: 
Do we know why? Proceedings of the National Acad- 
emy of Science, 99(10), 6524-6526. 

Meierding, T.C. (1982) Late Pleistocene glacier equilib- 
rium-line in the Colorado Front Range: A comparison 
of methods. Quaternary Research, 18, 289-310. 

Mishchenko, M.M., Dlugach, J.M., Yanovitskij, E.G., 
and Zakharova, E.T. (1999) Bidirectional reflectance 
of flat, optically thick particulate layers: An efficient 
radiative transfer solution and applications to snow 
and soil surfaces. Journal of Quantitative Sectroscopy 
and Radiative Transfer, 63, 409M32. 

Molnar, P., and England, P. (1990) Late Cenozoic uplift 
of mountain ranges and global climate change: 
Chicken or egg? Nature, 346, 29-34. 

Nakawo, M., and Young, G.J. (1981) Field experiments 
to determine the effect of a debris layer on ablation of 
glacier ice. Annals of Glaciology, 2, 85-91. 

Nakawo, M., and Rana, B. (1999) Estimate of ablation 
rate of glacier ice under a supraglacial debris layer. 
Geografiska Annaler, 81A, 695-701. 

Nakawo, M., Morohoshi, T., and Uehara, S. (1993) 
Satellite data utilization for estimating ablation of deb- 
ris covered glaciers. Paper No. 218 presented at Snow 
and Glacier Hydrology: Proceedings of the Kathmandu 
Symposium, IAHS Press, Wallingford, U.K., pp. 75- 
83. 

Nesje, A. (1992) Topographical effects on the equilibrium 
line altitude on glaciers. Geojournal, 27, 383-391. 

Nienow, P.W., Hubbard, A.L., Hubbard, B.P., Chandler, 
D.M., Mair, D.W.F., Sharp, M.J., and Willis, I.C. 
(2005) Hydrological controls on diurnal ice flow 
variability in valley glaciers. Journal of Geophysical 
Research, 110(F4), F04002.l-F04002.ll. 

Nolin, A.W., and Liang, S. (2000) Progress in bi-direc- 
tional reflectance modeling and applications for sur- 
face particulate media: Snow and soils. Remote Sensing 
Reviews, 18, 307-342. 

Nolin, A.W., and Stroeve, J.C. (1997) The changing 
albedo of the Greenland Ice Sheet: Implications for 
climate change. Annals of Glaciology, 25, 51-57. 

Nuth, C., and Kaab, A. (2011) Co-registration and bias 
corrections of satellite elevation data sets for quantify- 
ing glacier thickness change. The Cryosphere, 5(1), 
271-290. 

Osmaston, H. (2005) Estimates of glacier equilibrium line 
altitudes by the Area x Altitude, the Area x Altitude 
Balance Ratio and the Area x Altitude Balance Index 
methods and their validation. Quaternary Inter- 
national, 138/139, 22-31. 

Paul, F., Kaab, A., Maisch, M., Kellenberger, T., and 
Haeberli, W. (2002) The new remote sensing derived 


52 Theoretical foundations of remote sensing for glacier assessment and mapping 


Swiss glacier inventory: I. Methods. Annals of Glaciol- 
ogy, 34, 355-361 (Fourth International Symposium on 
Remote Sensing in Glaciology). 

Paul, F., Huggel, C., and Kaab, A. (2004) Combining 
satellite multispectral image data and a digital eleva- 
tion model for mapping debris-covered glaciers. 
Remote Sensing of Environment, 89, 510-518. 

Petzold, T.J. (1972) Volume Scattering Functions for 
Selected Ocean Waters (Tech. Rep. SIO Ref 71-78), 
Scripps Institution of Oceanography, San Diego, CA 

Prieur, L., and Sathyendranath, S. (1981) An optical 
classification of coastal and oceanic waters based on 
the specific absorption of phyplankton pigments, dis- 
solved organic matter, and other particulate materials. 
Limnology and Oceanography, 26, 671-689. 

Proy, C., Tanre, D., and Deschamps, P.Y. (1989) Evalu- 
ation of topographic effects in remotely sensed data. 
Remote Sensing of the Environment, 30, 21-32. 

Ramsay, B. (1998) The interactive multisensor snow and 
ice mapping system. Hydrological Processes, 12, 1537- 
1546. 

Raup, B., Kaab, A., Kargel, J.S., Bishop, M.P., Hamil- 
ton, G., Lee, E.M., Paul, F., Rau, F., Soltesz, D., 
Khalsa, S.J. et al. (2007) Remote sensing and GIS 
technology in the Global Land Ice Measurements from 
Space (GLIMS) Project. Computers and Geosciences, 
33, 104-125. 

Rees, W.G. (2006) Remote Sensing of Snow and Ice, 
Taylor & Francis, Boca Raton, FL. 

Robinson, D.A., and Kukla, G. (1985) Maximum surface 
albedo of seasonally snow-covered lands in the North- 
ern Hemisphere. Journal of Climate and Applied 
Meteorology, 24, 402 — 4 1 1 . 

Rossi, R.E., Dungan, J.L., and Beck, L.R. (1994) Kriging 
in the shadows: Geostatistical interpolation for remote 
sensing. Remote Sensing of the Environment, 49, 32-40. 

Ruddiman, W. F. (Ed.) (1997) Tectonic Uplift and 
Climate Change. Plenum Press, New York. 

Scambos, T.A., Dutkiewicz, M.J., Wilson, J.C., and 
Bindschadler, R.A. (1992) Application of image 
cross-correlation to the measurement of glacier 
velocity using satellite image data. Remote Sensing of 
Environment, 42(3), 177-186. 

Schaaf, C.B., and 21 others (2002) First operational 
BRDF, albedo nadir reflectance products from 
MODIS. Remote Sensing of Environment, 83(1/2), 
135-148. 

Seltzer, G.O. (1993) Late-Quaternary glaciation as a 
proxy for climate change in the central Andes. Moun- 
tain Research and Development, 13, 129-138. 

Shroder, J.F., Jr., and Bishop, M.P. (2000) Unroofing of 
the Nanga Parbat Himalaya. In: M.A. Khan, P.J. 


Treloar, M.P. Searle, and M.Q. Jan (Eds.), Tectonics 
of the Nanga Parbat Syntaxis and the Western Hima- 
laya (No. 170 in Special Publication), Geological 
Society, London, pp. 163-179. 

Siewert, C.E. (2000) A concise and accurate solution to 
Chandrasekhars basic problem in radiative transfer. 
Journal of Quantitative Sectroscopy and Radiative 
Transfer, 64. 109-130. 

Stroeve, J., and Nolin, A. (2002) New methods to infer 
snow albedo from the MISR instrument with applica- 
tions to the Greenland Ice Sheet. IEEE Transactions on 
Geoscience and Remote Sensing, 40(7), 1616-1625. 

Stroeve, J., Nolin, A., and Steffen, K. (1997) Comparison 
of AVHRR-derived and in situ surface albedo over the 
Greenland Ice Sheet. Remote Sensing of Environment, 
62, 262-276. 

Stroeve, J., Box, J., and Haran, T. (2006) Evaluation of 
MODIS (MOD10A1) daily snow albedo product over 
the Greenland Ice Sheet. Remote Sensing of Environ- 
ment, 105, 155-171. 

Torsnes, I., Rye, N., and Nesje, A. (1993) Modern and 
Little Ice Age equilibrium line altitudes on outlet valley 
glaciers from Jostedalsbreen, Western Norway: An 
evaluation of different approaches to their calculation. 
Arctic and Alpine Research, 25 , 106-116. 

Warren, S.G. (1982) Optical properties of snow. Reviews 
of Geophysics and Space Physics, 20(1), 67-89. 

Warren, S.G., Brandt, R.E., and O'Rawe, H.P. (1998) 
Effect of surface roughness on bidirectional reflectance 
of Antarctic snow. Journal of Geophysical Research, 
103(E1 1), 25789-25807. 

Williams, R.S., Hall, D.K., and Benson, C.S. (1991) 
Analysis of glacier facies using satellite techniques. 
Journal of Glaciology, 37(125), 120-128. 

Winther, J.G. (1993) Landsat TM derived and in-situ 
summer reflectance of glaciers in Svalbard. Polar 
Research, 12, 37-55. 

Zeitler, P.K., Meltzer, A.S., Koons, P.O., Craw, D., 
Hallet, B„ Chamberlain, C.P., Kidd, W.S., Park, 
S.K., Seeber, L., Bishop, M.P. et al. (2001) Erosion, 
Himalayan geodynamics, and the geomorphology of 
metamorphism. GSA Today, II, 4-8. 

Zemp, M., Haeberli, W., Hoelzle, M., and Paul, F. (2006) 
Alpine glaciers to disappear within decades? Geophys- 
ical Research Letters, 33, L13504, doi: 10.1029/ 
2006GL026319. 

Zwally, H.J., Abdalati, W., Herring, T., Larson, K., 
Saba, J., and Steffen, K. (2002) Surface melt-induced 
acceleration of Greenland ice-sheet flow. Science, 
297(5579), 218-222. 


